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ABSTRACT 


This thesis presents a study of the second-order effects 
in reinforced concrete frames. 

The load-deformation response of reinforced concrete 
Sway frames was obtained using realistic moment - axial load - 
curvature diagrams and the concepts of column deflection curves 
and sway subassemblages. The results of this analysis were used 
to study the App Bice beady of the moment magnifier and PA - methods 
of analysis. 

The parameters which have been studied include deflect- 
ions at ultimate load and service load, second-order moments, 
stability effects and the effective stiffnesses of reinforced 
concrete columns and beams. 

Both the moment magnifier method and PA - method were 
found to be applicable to material failures, but neither is satis- 
factory when dealing with stability failures. Approximate 
methods have been established to predict the failure mode. 

Deflections were found to exceed current code pres- 
cribed values in most cases. Stability failures exhibited the 


largest deflections. 
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Initial angle of CDC measured relative to the 
direction of the force F. 

Angle between tangent and the direction of the 
force F at the end of the j'th segment. 

Ratio of minimum end moment to maximum end moment. 
Positive for double curvature. 

Length of the j'th segment of the CDC curve. 

One half of story deflection. Inflection point 
assumed at mid-height of column. 

Deflection of the n'th floor relative to the floor 
below. 

Angle between the lines of action of the forces 

F and P. 

Shear distribution factor. 

Mom’ My 

Moy/M 

Ratio of reinforcement area to the area of the cross- 
section. 

Curvature of the j'th segment of the CDC curve. 


Mean curvature of the j'th segment. 
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CHAPTER I 
INTRODUCTION 


In recent years the use of slender buildings and 
building components has become very common, and from a structural 
point of view this trend has made it necessary to consider certain 
aspects of structural behaviour that are not important when 
members and structures have a low height to width ratio. 

Of particular interest in the analysis of building 
frames are the deflections produced by the lateral loads. Phe 
deflections can be relatively large when the building has a large 
height to width ratio. The axial loads acting in the columns will 
produce additional moments and forces in the structure when acting 
through the lateral displacement of the columns. When these 
effects are taken into account in the structural analysis of the 
building, the analysis is referred to as a "second-order analysis." 
Fig. 1.1 (a) shows a simple frame acted upon by a lateral load Q 
and vertical loads P applied at the top of the columns. The de- 
formed shape of the frame is shown in Fig. 1.1 (b). From this 
figure it is seen that when the structure has deflected laterally 
the vertical load P will also contribute to the lateral sway of the 
frame. When the frame has reached its final deflected position the 
axial load produces a sway moment PxA which is commonly referred 
to as the "PA-moment". 

The purpose of this thesis is to study the effect of the 
PA-moments on the Peeve of reinforced concrete frames. A compu- 


ter analysis is used to generate the load-deformation behaviour of the 
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SECOND-ORDER EFFECTS ON SWAY PERMITTED FRAMES 
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analytical model and the results of this analysis are used to study 
the applicability of current design methods and a proposed method 
of second-order analysis. 

The general theory which forms the basis of the computer 
analysis is presented in Chapter III. Moment-axial load-curvature 
curves are generated to establish the load-deformation behaviour 
of the member cross-section and the concepts of column deflection 
curves and sway subassemblages are also presented. 

Chapter IV presents the analytical model used in the 
computer analysis and load-deformation curves of the model. 

These curves were computed using the concepts presented in Chapter 
III. The significance of the deflection indices obtained in the 
analysis is discussed and effective elastic column stiffnesses are 
presented. Finally, the use of the sway subassemblage concept in 
frame analysis is discussed. An extensive set of sway subassemblage 
charts are presented in Appendix A. 

Effective elastic EI values for reinforced concrete beams 
are studied in Chapter V and recommended design values of flexural 
stiffness are presented. 

Chapter VI is devoted to the discussion of stability 
problems in multi-story frames. Results from current methods of 
second-order analysis are compared to the results obtained in the 
computer analysis. An alternate method of analysis is also pro- 


posed and its applicability and limitations are discussed. 
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‘CHAPTER II 
REVIEW OF CURRENT METHODS OF SECOND-ORDER ANALYSIS OF BUILDING FRAMES 


2.1 Introduction 

The importance of a realistic second-order analysis has 
become evident in recent years due to the use of smaller load 
factors and more slender compression members in building frames. 
Deflection limits based on serviceability requirements have 
become an important design criterion in current building codes and 
inmany cases are a governing factor in the design of tall buildings. 

Although large deflections frequently do not affect the 
strength of beams or low structures, this is not so in frames which 
are free to sway under lateral loads. The lateral deflections give 
rise to additional moments and may produce a significant reduction 
in the load carrying capacity of the structure. It is, therefore, 
essential that these second-order effects can be computed in order 
to establish a rational design procedure, and considerable research 
has been carried out in recent years to improve existing methods and 
to develop new methods of accounting for second-order effects. 

Due to the complexity in carrying out an exact second- 
order analysis, approximate methods of various types have been dey- 
eloped for design office use. Most building codes suggest design 
methods for compression members which are of a semi-empirical 


nature or are based on approximate theorectical solutions. 


oe Approaches to Frame Stability Problems in Building Codes 


There are two main methods of considering the stability 


of concrete columns and frames in use to-day. One is the moment 
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magnifier method presented in ACI 318-71 (1,2,3) 


and the other is 
the complementary moment method adopted by the CEB (Comité Européen 
du Beton) (425). 

In the moment magnifier method the moments obtained from 
a first-order analysis are multiplied by a magnification factor 
which is a function of the ratio of the applied axial load to the 
critical load of the column. To account for the end restraints on 
columns in rigid frames an effective length factor k is introduced. 
The restrained column is then replaced by a pin-ended column of 
length kl which is designed for the applied axial load andthe 
magnified moment. Nomographs (2) obtained from an elastic solution 
may be used to determine k for columns which are assumed to be 
either free to sway or fully braced. The accuracy of the method 
depends to a large extent on the accuracy of the relative values of 
column to beam stiffness. The ACI Code provides expressions to 
compute the column stiffness as a function of the stiffness of the 
gross concrete section EQlg: and the stiffness of the reinforcement 
El. The variation in stiffness due to the axial load is not in- 
cluded. 

The CEB method, which has been adopted in one form or 
another by many European countries, accounts for second-order 
effects by adding a complementary moment to the first order moment. 
The CEB Recommendations present a simplified method for computing 
the complementary moment which is given as a function of the effect- 
ive length and approximate expressions for the curvature. The 


effect of rotational end restraint on the column is accounted for in 
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a manner similar to that in the moment magnifier method. The 
actual column is replaced by an equivalent pin-ended column of 
length kl which is then used to determine the complementary moment. 
The effective length factor is assumed to be obtained from a linear 
elastic analysis and is therefore subject to the same assumptions 
as used in establishing beam and column stiffnesses. No specific 
method is suggested in the CEB recommendations for evaluating k. 
The calculation of effective lengths is discussed more 


fully in Section 6.3 of this thesis. 


2.3 Iterative Procedures 
Several authors have suggested the use of iterative pro- 
cedures, often referred to as P-A methods, to carry out second- 


(6575859510) the basic idea behind these 


order frame analysis 
procedures is that the moments produced by the total vertical load 
P, acting through a lateral deflection A at a certain level ina 
building may be replaced by equivalent lateral shears applied at 
floor levels. 

The equivalent shears are computed from the deflections 
obtained in a first order analysis and added to the applied 
lateral loads. This process is repeated until the deflections in one 
cycle agree, within specified limits, with the deflections obtained 
in the previous cycle. 

A rapid convergence of the process is essential to make 
it feasible for design office use, thus no more than three iterations 


should be necessary. When the axial load has a value close to the 


critical load the convergence tends to become very slow. 
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K. Aas-Jakobsen has proposed a finite element approach 

to solve for secon-order effects under linear-elastic conditions. The 
stiffness matrix [K] is assumed to be the sum of two stiffness matrices 
[Ky] and LK» J where [Ky] is the first order stiffness matrix and [Kp] 
is the non-linear second-order stiffness matrix. The matrix [Ko] is 
obtained through an iteration procedure. When unit displacements are 
applied to the member the axial load required to maintain equilibrium is 
unknown and can only be obtained by trial and error. Aas-Jakobsen 
suggests that the axial load be set equal to zero in the first cycle. 
From the first order forces obtained in the first cycle an equivalent 
axial load can be computed and used in the second cycle. The process 

is repeated until the axial load found in one cycle is close 
to the value computed in the previous cycle. Aas-Jakobsen 
states that the process will usually converge rapidly so 


that two cycles are generally sufficient. 


2.4 Non - Linear Analysis 


The majority of today's structures are designed 
using a straight line force - deformation relationship. The 
validity of such an assumption depends largely on the degree to 
which non-linearities affect the distribution of forces. The 
moment in continous beams are relatively insensitive to the 
value of EI used unless the relative EI changes during the 
loading history. Building codes usually specify EI values that 
will give results on the conservative side. On the other hand, 
in second-order analysis where the actual lateral deflections 


are important, a reasonably good approximation to the actual 
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FIG. 2.1. SIMPLIFIED NON-LINEAR M-¢ DIAGRAMS 
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load-deformation curve is necessary to properly predict the 
effect of lateral deflections. 

Structural analysis based on a load-deformation curve 
which is approximated by a non-linear smooth curve or by several 
straight lines are usually referred to as "non-linear analysis". 

Due to the difficulty in representing the load- 
deformation curve of reinforced concrete cross-sections by an 
analytical expression, straight line approximations have been 


(12) 


suggested It has been found that three straight lines are 


sufficient to give results close to test results (13,14,15) 

Fig. 2.1 shows some typical straight line aoiuen baer vetane dia- 
grams for various members and loading conditions. The points 
C,YT,YC and F represent cracking of concrete in tension, yielding 
of the tension reinforcement, yielding of the compression rein- 
forcement and failure due to crushing of the concrete, respectively. 

The region o-c represents the uncracked state. When 
cracking has occurred, region C-YT, there is a loss of stiffness 
but the member is still behaving elastically. There is a further 
reduction in stiffness after yielding of the reinforcement in 
YT-F, and the member no longer behaves elastically. 

A non-linear analysis of complex structures requires 
considerable computation and its use will usually not be justified 
in practice. Hence, it is desirable, in the case of reinforced 
concrete structures, to find elastic stiffnesses that lead to 
acceptable results or to find other ways of carrying out the 
structural analysis. Both techniques are explored more fully 


later in this thesis. Section 4.5 in Chapter IV presents an 
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investigation into effective EI values for reinforced concrete 
columns and in Chapter V the same problem is considered for rein- 
forced concrete beams. Chapter IV also presents the development 
of an alternate method of analysis as applied to reinforced con- 
crete structures. 


Breen (16) 


used a numerical integration process to 
compute the relationship between end moment and end rotation for 
reinforced concrete members taking into account the non-linearity 
of the moment-axial load-curvature diagram. By plotting the 
individual beam and column moment-rotation curves on the same 
graph he showed that the moment-rotation curve of a frame joint 


could be established and from these curves the loading capacity 


of the frame was predicted. 


2.5 The Sway Subassemblage Concept 

A semi-graphical method for designing girders and 
beam columns in rigidly jointed unbraced frames has been proposed 
by Daniels and Lu (17), 

The procedure starts by analyzing a frame whose members 
have been selected in a preliminary analysis. Each story is 
subdivided into a number of sway subassemblages consisting of a 
column and one or two adjacent beams, depending on whether the 
column under consideration is an exterior or interior column. 

To trace the load-deformation behavior of each sub- 
assemblage a set of restrained columns design charts is used. 
These are non-dimensional plots of lateral load against lateral 


deflection of half a column for various values of end restraint, 


slenderness ratio: and axial load and include the effect of 
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lateral deflections on the load carrying capacity of the column. 
When the load-deflection curve has been established 
for each of the subassemblages in a story for a given axial load 
level, the curves are combined to give the overall load-deflection 
curve for that story for that axial load level. The lateral 
load at which the slope of the curve for the story becomes 
negative corresponds to instability of the story. In this 
manner the load-deflection relationships of a story may be deter- 
mined without prior design of other parts of the frame. 
In Chapter III the development of sway subassemblage 
charts for reinforced concrete columns will be discussed in more 


detail. 
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CHAPTER III 


DEVELOPMENT OF COLUMN DEFLECTION CURVES AND SWAY SUBASSEMBLAGE 
CHARTS FOR REINFORCED CONCRETE COLUMNS 


3.1 Introduction 

In the analysis of a building for lateral loads each 
story of the frame can be divided into a series of subassemblages 
each consisting of a half column and adjoining beams. The sub- 
assemblages can be analyzed separately and their individual 
lateral load-deflection curves superimposed to obtain an 
overall lateral load-deflection curve for the story. Sway 
subassemblage charts consisting of a series of restrained column 
curves for various degrees of end restraint have been developed 
to aid in the analysis of structural steel subassemblages as 
described in Section 2.6. The development of such charts for 
reinforced concrete columns will be discussed in this chapter. 

The properties and use of these charts and curves will be dis- 
cussed in Chapter IV. 

The sway subassemblage charts and restrained column 
curves are developed using Column Deflection Curves, usually 
referred to as CDC curves. A CDC curve is a plot of the deflected 
shape of a column subjected to a given axial load with a given 
slope at the point where the moment equals zero. Such a curve may 
be used to represent the deflected shape of a beam-column having the 
same axial load, but different end conditions, which may vary from 
hinged to fully fixed. If a family of such curves exists 
for a certain axial load it is possible to determine the 


relationship between the applied lateral load and the corres- 
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ponding lateral deflection that must exist for the system to 
be in equilibrium under given end restraints. 

In order to compute a CDC curve it is necessary to 

Know the response of the column cross section to external 
loads, which, in this case, is conveniently represented by 
Moment-Axial Load-Curvature curves (hereafter referred to as 
M-P-¢ curves). The shape of these curves is a function of 
material properties, geometry and external loads. 

The steps involved in establishing the sway subassem- 

blage charts may be summarized as follows: 

1. Determine material stress-strain relationships, 

2. Compute the M-P-¢ curve for the cross-section 
for a given axial load, 

3. Generate a family of CDC curves for the given 
axial load, 

4. Compute the sway subassemblage charts for the 
given axial load and various values of end 
restraint. 

Each of these steps will be discussed more fully in the 


following sections. 


3.2 Material Properties 
3.2.1 Stress-Strain Curves for Concrete 

The stress-strain relationships assumed for concrete in 
compression and tension are shown in Fig. 3.1(a). The compression 


curve is similar to the one proposed by Hognestad (18) 
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FIG. 3.1(a) | STRESS-STRAIN DIAGRAM FOR CONCRETE IN FLEXURE 
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FIG. 3.1(b) | REINFORCEMENT STRESS-STRAIN DIAGRAM 
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The non-linear part of the compression curve is given 


by Equation (3.1) and the stress-strain relationship for concrete 


in tension is represented by Equation (3.2) (19), 





where 


The value of the modulus of elasticity of concrete has 


been taken as 


5 


E, = 18 x 10° + 500 f, psi 


The compressive strength of concrete loaded in flexure, 


un was assumed to be 85 percent of the strength of control Bae 


and under the same loading conditions the ultimate strain was taken 
as 0.004. 
The split cylinder strength fi, is assumed to be equal 
to 7.0 f° and the ultimate tensile strain, Euier May be expressed as 
1 
2t/E.. 


3.2.2 Reinforcing Steel Stress-Strain Curve 


Only reinforcing steel with fy = 60000 psi has been 


considered in this study and the assumed stress-strain curve for this 
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study and the assumed stress-strain curve for this steel is 

shown in Fig. 3.1(b). It consists of three parts; the initial 
elastic region, a flat yield plateau and a strain hardening 
region. The length of the flat plateau has been assumed slightly 
longer than what is usually indicated by tests. The effect of 
this is to overestimate somewhat the deflections at ultimate. 
This in turn results in an overestimate of the second order 
moments and therefore tends to be conservative. The slope of 

the strain hardening portion was based on the average value from 


a number of available stress-strain curves. 


3.3 Moment-Axial Load-Curvature Relationships 
3.3.1 Calculation of M-P-o Curves 


The computation of the M-P-¢ curves is based on the 
following assumptions: 
a) Plane sections remain plane. 
b) The stress-strain diagram of concrete 
is as described in section 3.2.1. 
c) The concrete is assumed to crush when the 
strain in the extreme compression 
fibre reaches 0.004. 
d) The concrete cracks when the tensile strain 
exceeds the ultimate tensile strain. 
e) The stress-strain diagram for steel is as 
described in Section 3.2.2 and shown in Fig. 3.1 (b). 
For a given axial load the moment-curvature relationship can be 


computed using a trial and error procedure to find the neutral axis. 
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This process is rather tedious for hand calculations and a computer 
program (20) was therefore used. 

The program can handle any shape of section provided 
there is bending about one axis only, and includes the effect of 
tensile stresses in the concrete and strain hardening in the rein- 
forcement. The basic steps in the computations are as follows: 

1. Divide the cross section into a number of elements 
of sufficiently small size such that the variations 
in strain across the element is relatively small. 

2. Assume a small value of the curvature and a position 
of the neutral axis. 

3. Compute strains, stresses and forces for each element. 

4. Check if the forces acting on the cross section are in 
equilibrium. If they are not, choose a new position of 
the neutral axis and repeat. 

5. When equilibrium has been established, compute the 
internal moment. 

6. Repeat 1 through 5 for increasing values of curvature 
until the ultimate concrete strain in compression has 


been reached. 


3.3.2. The Effect of Various Parameters on the M-P-¢ Curve 


The major variables affecting the shape of the M-P-¢ curves 
include: 
a. The axial load ratio Pu/Po (or Pu/P, ) 


b. The steel percentage, .p 


c. Material strengths, fa fy. 
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d. Cover ratio, y 

e. Shape and size of cross-section 

f. Distribution of reinforcement over the 

cross-section. 

Only axial load, shape of cross-section and distribution 
of reinforcement will be considered here. 

Fig. 3.2 shows four curves computed for the same columns 
with four different values of axial load. 

As the axial load is increased the failure mode changes 
from ductile failure, shown by point 4 in this figure, to brittle 
failure ,point 11. The increase in axial load also delays the opening 
of cracks on the tension side of the column. Initiation of cracking 
occurred at points 1,5 and just above 8. The heavily loaded columns 
experience little reduction in stiffness until the compression steel 
yields at 8 and 10. Yielding of the tension steel occurs at 2 
while 6 represents yielding of the compression steel followed by 
cracking of the concrete. Points 3,7,9 and 11 represent points of 
maximum moment capacity for P/P. equal to 0.2, 0.4, 0.6 and 0.8 
respectively. 

In Fig. 3.3 the M-P-¢ relationship has been plotted for 
the three different types of column for the same axial load ratio. 
The cross-sectional properties are given in Figure 4.1. All three 
columns had the same total steel percentage and the same cover. 

All sections show practically the same amount of rotation capacity 


at ultimate but the moment capacities are reduced considerably for 
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type 2 and type 3 due to the less efficient placing of the 
reinforcement. Some of the reduction for type 3 is also due to 
less efficient distribution of the concrete. 

Distributing the reinforcement over the cross-section is seen 
(Fig. 3.3. ) to produce a much smoother transition from the elastic 


to the inelastic range. 


3.4 Column Deflection Curves 
3.4.1 Introduction 


In 1910 von Karman used CDC's to determine the strength 


(21), 


of beam-columns with small eccentricities Chwalla generalized 


von Karman's work in 1934 to include other loading conditions (22) | 


Since then several investigators have expanded Chwalla's basic 
ideas in various directions (23 to 29) | 
The basis of the use of CDC's as applied in this thesis 


was established by Ojalvo (27), 


3.4.2 Basic Relationships of Column Deflection Curves 

A typical beam-column under applied loads is shown in 
Fig. 3.4(a). The equilibrium equations of this member may be 
established in the usual manner from statics, and the external 
moment at any cross section a distance x from the left hand end is: 

Me = M jist (s-1)]+ Pv. (3.3) 
where Vi is the deflection at distance x and g is the end moment 
ratio, positive for double curvature. 

The external moment must equal the internal moment. Thus: 


M. = M, = f(g) = f(-v,) (3.4) 
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FIG. 3.4 FORCES ON A BEAM-COLUMN 
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or: 

M E +7 (8-1) | + PV = f(-v,) (3.5) 
where f(¢), etc. means a function of 9, etc. 

However, a different approach may be used to express 
equilibrium. The system of forces in Fig. 3.4(a) may be replaced by 
a single force F as shown in Fig. 3.4(b). 

The direction of the force is defined by the angle y such 
that 

Fa = M and Fb = eM, 
and in general 


M = Fv = f(-v ) 


E (3.6) 
Summation of moments about A and B taking counterclockwise 
moments as positive yields, respectively, 
Fl(m+1) sin (-y) = -M (3.7) 
and 
Fmi[sin (-y)]= eM (3.8) 
Combining Eqns. (3.7) and (3.8) gives the point of intersection 


between the line of action of F and the line AB, 


2 a6 
bie BE (3.9) 





To satisfy the axial load requirement it is necessary that 


joc (3.10) 
COSY 


and by combining Eqns. (3.8),(3.9) and (3.10) the angle y may be 


determined from 


tan y = (1+8)M (3.11) 
Pl 
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Thus Eqns. (3.9),(3.10) and (3.11) define completely the 
magnitude, location and direction of F, and the equilibrium equation 
is now Equation (3.6). 

To obtain the deflected shape of the beam-column it is 
necessary to integrate Equation (3.6), and the M-P- curves form the 
basis for this integration. The process does not have to be confined 
to the section AB, but may be carried beyond these points. This is 
illustrated in Fig. 3.5 where the curve has been extended to intersect 
the line of action of F. The beam-column AB is a part of this half 


wavelength whose length and initial angle are L and Os respect- 


CDC 
ively. Just as the portion AB of the CDC represents the equilibrium 
of the beam-column under discussion, so another portion of the CDC 

is the shape of another beam-column. Thus any one CDC can give 
information about an infinite number of beam-columns with various 

end conditions. Fig. 3.6 shows some examples of how the CDC may be 
used to represent various types of beam-columns. Column AB is hinged 
at A and has applied moment at B. CD is bent in double curvature with 


end moments of equal magnitude and EF is a column bent in single 


curvature with unequal end moments. 


3.4.3 Assumptions in the CDC Calculations 


Since the M-P-# curves form the basis for CDC calculations 
it follows that the CDC's are subject to the same assumptions which 
were used in developing the load-deformation response of the cross 
section. 

Furthermore, it is common to assume that the equivalent 


axial force F is equal to the actual axial load P,. From Eqn. (3.10) 


it is clear that this assumption holds for small values of y. However, 
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FIG. 3.6 CDC'S REPRESENTING VARIOUS TYPES OF BEAM-COLUMNS 
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Eqn. (3.11) indicates that for a beam-column of length 1 and tid 
moments M and gM, tan y may become large for small values of axial 
load and non-zero values of (1 + 8). If the end moments are equal 
and opposite (g= - 1, single curvature) y = 0, F = P and the solution 
is always exact. But as the ratio of end moments increases to a max- 
imum of + 1, i.e. perfect double curvature, tan y increases to a 
maximum of 2M/P1 The worst situation arises for a short column 
with a low axial load and a maximum moment equal to the ultimate 
moment. Consider the values from Table 4.1 for a load ratio of 
0.2 and a slenderness ratio 1/h=70. From Eqns. (3.11) and (3410) 
it is found that 

F = 1.04P 
Thus by using F = 0.2P in the analysis the column has in fact 
been analyzed for an axial load which is less than the actual value. 
In this case the discrepancy is only about 4% but for shorter 
columns and for smaller axial load the error would be more severe. 
For practical purposes, however, this is not likely to be a serious 
problem since most columns will not have a combination of P/P, and 
1/h within the critical range. 
3.4.4 Calculation of CDC Curves 

The basic principles underlying the concept of column 
deflection curves were outlined in Section 3.4.2. In this section the 
equations required to compute the shape of the CDC curve will be pre- 
sented. 

To account for the distribution of cracking along the 
member it is necessary to divide the curve into a number of segments. 


The length. of the segment should be such that the degree of cracking 
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(30) 


is nearly constant within the segment. Galambos has suggested 
that a segment length equal to 4r will give good results for 
steel members. A somewhat shorter length would probably be appro- 
priate for reinforced concrete members since 4r generally exceeds 
the thickness of the column. In the calculations reported in this 
thesis the segment length was arbitrarily chosen as 5 inches. 

A typical segment of a CDC curve is shown in Fig. 3.7. 
If the segment is short it may be assumed that the curve is a 
circular arc. The radius of the j'th segment is then 1/$;. 

From the geometry of Fig. 3.7 the following relationships 
may be established (ignoring second-order terms such as aby. 


2 Be 
a4 and oy): 


bd = a 54 
54¢ 
angle dac = a 
ac = 6, 


ou2, 
cd 879,42 
The deflection and the angle of the tangent at the end of the first 


segment may now be written as 


2 
a6) 7 (674,/2) (3 12) 
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ST ee 649) (3.13) 


And in general at the end of the j'th segment 
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a, = a, - 6.9 


oe Ys aes BR (3.15) 
The values of Vs and a; are now used to determine V54y and A547" 
The first step in the calculations is to determine a 
value of O and 6. The length need not be the same for all segments 
but usually is for convenience. The value of $ may be determined 
from the M-P-¢ diagram when the moment has been determined. This 


moment is taken as the mean moment in the segment and the general 


expression for the mean moment in the j'th segment is 


Maj 7 Ysa * F 8j%5-3 (3.16) 

By entering the M-P-o curve with the value of a the 
curvature of the j'th segment may be determined. 

Due to the difficulty in obtaining a relatively simple 
analytical expression for the M-P-¢ diagram, this curve was approxi- 
mated by a number of straight lines. A high degree of accuracy may 
be obtained by using a sufficient number of lines. A typical example 
of such an approximation is shown in Fig. 3.8. 

If it is assumed that M. < Ming < Mega where the subscript i refers 
to the i'th point on the M-P-¢ curve, then the curvature corres- 
ponding to Maj is given 


by 
M 


.-M 
ba mJ 
omnj ey? t 


F 
6.475 facds) 
mad | tI \ 


The value of $nj is now substituted for F into Eqns. (3.14) and 





(3.15) to compute the deflection and the slope of the tangent, res- 
pectively, at the end of the j'th segment. 
The computations may be terminated when one quarter of the 


curve has been computed since the curve is symmetrical. When Mag 
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FIG. 3.8 LINEARIZED MOMENT-CURVATURE DIAGRAM 
- TIED COLUMN WITH BARS IN FOUR FACES, Pio = 0.4 
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reaches the moment capacity of the cross section the above equations 
are no longer valid because they do not account for angular dis- 
continuity at the hinge. 


3.4.5 Comparison of CDC Calculations with Test Results 


The theoretical calculations have been compared to tests 


carried out by Breen (6) (31) | 


and Chang The CDC's corresponding 
to failure of the columns are plotted with the test results in 
Fig. 3.9 and Fig. 3.10. 

The test specimens had heavy end brackets and hence the 
section properties were not uniform along the length of the member. 
To get a meaningful comparison it was necessary to ignore the end 
brackets and measure deflections from the line joining the ends of 
the actual column. This part of the member is represented on the 
CDC by the distance AB. The test results were plotted using AB as 
the horizontal axis allowing for the deflection at the end of the 
bracket. 

The maximum deflections from the CDC calculations are 
about four percent below Breen's test result and about ten percent 
above Chang's result. Thus the deflected shape of these columns 


was predicted with reasonable accuracy from the column deflection 


curve. 
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3.5 Sway Subassemblage Charts 
3.5.1 Basic Equations 

The basic ideas behind the subassemblage concept were intro- 
duced in Section 2.6, and the tools required to apply the concept, 
i.e. M-P-o diagrams and CDC curves, were presented in Sections 3.3 
and 3.4, respectively. It is clear from those discussions that 
considerable computational work is required if this method is used 
from first principles and that it would be quite unsuitable for 
hand calculations. A more practical approach is to develop stan- 
dard design curves which are applicable to a large number of cases. 
Such curves have been developed to aid in the design of steel 
structures(27), 

A multi-story frame is shown in Fig, 3.11(a). The sub- 
assemblage consisting of the n'th floor girders between column lines B 
and D and half of the column in line C between floors n and n+l is shown 
in Fig. 3.11(b). A point of contraflexure has been assumed at mid- 
height of the column. 


The model is acted upon by a lateral load 


Qh =A H. (3.18) 
Where x is a column shear distribution factor, the axial 

load Es which is the sum of all vertical loads acting on column line 

C above level n and the moment Ma-1 which is produced by the column 

shears above the n'th floor. A further simplification has been 

made in Fig. 3.11(c) where M. represents the total restraining 

action supplied by the girders. The angle e is the rotation of 


the joint and y is the angle between the chord and a tangent to the 
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upper end of the column curve. y is determined from the CDC curves 


and at the end of the j'th element it is given by equation (3.19). 
V. 
a 
ome F a, (3.19) 
The moment at the upper end of the column is given by the 
equilibrium Equation (3.20). 


4 A 
sya iba cae (3.20) 
Un a GOD antiynie2 





Equilibrium of moments at the upper end requires that 








M. = My + Mind (3.21) 
The moment My-1 may be expressed as 
] A 
mn n-1 n-1 (3.22) 
Dee ste i nate 
Since Qn-1 < Qn and Pel < C the value of M will always be 
smaller than My unless An-1 is appreciably greater than Ay 
Normally, ae ee and Ma-1 may be assumed to be equal to Mae 
This assumption will then be conservative. 
Equation (3.20) may be rewritten as 
M. = aM (3.23) 
The relationship between the rotations y, 6 and 4 a 
is found from the geometry of Fig. 3.11(c) and is given by 
An 
Taare sae ¥ (3.24) 
n 


If the restraining action provided by the beams is assumed to be 
linear the restraining moment may be expressed as 

M. = KeM (3.25) 
where K is the beam stiffness and yl is the ultimate moment 


_ capacity of the column. 
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sta. Procedure for Computing Subassemblage Curves 


The computations for a given subassemblage curve start 
by specifying a value of the axial load and the restraining function. 
For various values of 6 the restraining moment is calculated from 
Equation (3.25) and the moment at the top of the column is deter- 
mined from Equation (3.23). When the rotation y has been found, 
usually from a CDC analysis, the story deflection is determined 
from Equation (3.24). The lateral load which is compatible with 
the joint rotation and the forces acting on the column is found 
from Equation (3.20). 

The results can conveniently be plotted in a non-dimen- 
sional form of lateral load against story deflection. To accomp- 
lish this Equation (3.20) is written in the form given by Eqn. (3.26) 
and Qa! My is plotted against d/l) 


An 


aha 
a ene a od dey (3.26) 
2M, M 


A computer program was written to generate the restrained 
column curves and a description of this program is presented in 
Appendix B. The properties and use of these curves are discussed 


in Chapter IV. 
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- CHAPTER IV 
RESTRAINED COLUMN CURVES FOR REINFORCED CONCRETE COLUMNS 


4.1 Introduction 
. The load deformation response of a sway subassemblage 

can be described by means of restrained column curves which represent 
the load-deformation characteristics of a column having a given 
rotational end restraint. The load-deflection relationship of a 
subassemblage at any given stage of loading is obtained from the 
restrained column curve corresponding to the same end restraint as 
that provided by the beams of the subassemblage. Section 4.7 
describes in more detail how restrained column curves are used to 
obtain the complete load-deformation response of a sway subassemblage. 

The general theory relating to restrained column curves 
and sway subassemblages was presented in Section 3.5. In this chapter 
the theory will be used to generate restrained column curves for 
reinforced concrete columns. The series of curves, contained in 


Appendix A were generated by the computer program in Appendix B. 


4.2 The Analytical Model and Variables Studied 
To study the general properties of subassemblage charts 
for reinforced concrete members and to provide comparison with other 
design methods, charts were generated for three columns: 
Type 1. square column with reinforcement in two faces, 
Type 2. square column with reinforcement distributed in 
all four faces, 
Type 3. circular column with reinforcement distributed along 


the circumference. 


39 


; vs u@ ak 
» — 




















i Te : é 
’ 

i ~ ; is 
e 7 ie , 
Pa aoe 

v4, 
WGA HOY VAWRUD iMMlJOD GARTANTZIA 
ad - TV a* Me id \. ' : 7 
ay: 
i 
notisubotsel 1" 


mr gaToqany ‘ 
Pa PAties e 





at 


H 
- 


a 





y sel 
ae : -i ~~ 8 
+ | hs : a v : 
o i ber é 
=¢3 — ab 13 Bret toe et. ‘” t 7 


¢ ite 
PSN % SHASM 4G be if ‘928 od nés> ¢ 
t - ~ ’ < 


Ji terveyos verte shottsmetsb-bsol edd 7 


-" 


—_ > | 

bsof atT .anbettess bas [snore ssor 
~ 

f 7 


’ . - ae _ - 4 4 
st2 tievie yns 26 apeldmeczadue - 
‘ 


) ~ otsb-béc! 9 
2) 
0 f ' of ; ey vrouT lsyanst at 
"S7gen< ifnl  .c.&° af. if bsdIns2std> 2hw 266 fi Seoeadlais vane m8 
. J | Po 
2oyTUD nuuloo berntsyfes steyonsp oF- ‘beau aa Ffiw yeas odd 


“= ; bad 


itBinoo .2evwo. %o 29fte92e sfT .2him aati oo, ays" 0% beovotnts 










4 xtbpsagé mt meteor vecugmoo enti vd seseranse, ‘o18w A bneqs in 
: hae 


r. 


= 





buf? gsidstw bos ‘fobo [aottyl sna 


2dvedo opsldmezesdyu2 to eotirasqorg lasanap aris ybuse a A 


® 
i 

‘s 

2 





(ato ftiw qoztraqmes sbfvevg of bins jrodeam steran09 beototnte" 


= hic: 
‘anmifoo seit 107 boderonse care etusds ceboridom ng teab 


e 


a Se 
- 29087 owt nt sSnemsoyvetats? ride eater Some Ls ; aT 
- ft bedadindel 


The column cross-sections are shown in Fig. 4.1 along with the mater- 
ial properties. 

The subassemblage considered in this analysis is shown in 
Fig. 4.2(a). The springs at the ends of the beam and the column re- 
present the net restraint provided by the frame which the subassemblage 
is assumed to be a part of. It will be assumed that the stiffness of 
the springs is such that the members are bent in double curvature 
with a point of contraflexure in the middle. Thus the model may be 
represented by Fig. 4.2(b). 

It will be assumed that no upper column frames into the 
beam-column joint. Hence, Equation (3.23) may be modified to read 

M. = M (4.1) 

The equilibrium equation for the column may be written down 
from Fig. 4.2(c) 

M= QL + Pa, (4.2) 


and in non-dimensional form 


(4.3) 


—|c 
— 


QL M-P 
M 


i M 


Cc 


The lateral load vs. lateral deformation relationship may be 
traced by plotting QL/M,, against A/L. 

The model was analysed for five values of beam stiffness K 
(see Equation (3.25)): 100, 200, 400, 600, and infinity. The values of 
K derived in Section 5.3 for a representative T-beam and flat plates 
are about 500 - 2000 and 50 - 140, respectively, when combined with the 
columns shown in Fig. 4.1. | 

Eight values of the slenderness ratio 1/h were used, varying 


from 5 to 40 in steps of 5. 
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Pp (a) TYPE 1 
Ff. = 4000 psi 
y #60000 psi 
10" (b) TYPE 2 
p =0.02 
Y =O 
COVER 1.25 in 
10" cel YPE. 3 


FIG. 4.1 COLUMN CROSS-SECTIONS 
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FIG. 4.2 ANALYTICAL MODEL FOR DEVELOPMENT. OF 
RESTRAINED COLUMN CURVES | 
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The complete set of charts obtained from this analysis is presented 
in Appendix A. A total of 146 charts are presented for 690 indivi- 
dual combinations of length, shape and end restraint. 

The coordinates of the linearized M-P-¢ diagrams used 


in the analysis are given in Tables 4.1, 4.2 and 4.3. 


4.3 Basic Properties of Subassemblage Curves 


Some typical results from the analysis are shown in 
Figs. 4.3 and 4.4, and the moment curvature diagram used in the calcu- 
lations of these curves is shown in Fig. 4.5. | 

The line ab in Fig. 4.3 represents the locus of columns 
failing as a result of material failure, due to reaching the cross- 
sectional capacity at the top of the column. This line slopes be- 
cause the PA moments reduce the lateral load capacity. Lines parallel 
to ab give the capacity of the subassemblage when a plastic hinge 
forms in the beams prior to the column reaching its capacity. Thus 
the capacity of the subassemblage is limited by the line cd when 


the restraining capacity of the beams is 0.75 Mae 


The point A in Figs. 4.3 and 4.4 corresponds approximately to 


point 1 on the M-P-¢ curve shown in Fig. 3.2. This is the point 
where the column starts to crack and there is a loss of stiffness, 
but it continues to behave elastically. 

While the short column in Fig. 4.3 reaches its ultimate 
capacity with a corresponding maximum value of Q1/M, the longer 
column in Fig. 4.4 exhibits material failure under a lateral load 
which is less than its load carrying capacity. The practical limit 


of this column is reached on the line mn. For any increase in load 
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Axial load ratio, P /P 
fos [ous 
, 2 
0. 


F ie) 
226.5 1315.0 |360.0 | 410.0 | 410.0 0 | 210. 
¢x10° | 0.071 |0.100 |0.120 10.147 |0.158 3 10. 





u 
123 | 0.093 | 0.085 
6x10° 0.424 {0.310 | 0.136 | 0.163 | 0.193 | 0.197 | 0.162 | 0.114 
6x10° 0.790 | 0.416 | 0.218 | 0.186 | 0.270 | 0.246 | 0.197 | 0.153 
$x10° 1.190 | 0.491 | 0.295 | 0.245 | 0.329 | 0.269 | 0.288 | 0.182 
Moment | 671.0 | 749.0 | 663.0 | 663.0 eed bee 299.0 
6x10° 1.880 | 0.557 | 0.365 | 0.341 | 0.377 | 0.326 | 0.296 | 0.200 


ee ea 
4x10° 0.710 | 0.500 | 0.424 | 0.464 0.212 
Seer oe 
$x10° 0.876 | 0.525 | 0.533 

i ee 
$x10° 0.592 | 0.647 


Units: moment is inch-kips, curvature is 1/inch. 










TABLE 4.1 


Coordinates for linearized M-P-¢ diagrams for Type 1 column. 
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Axial load ratio, ae 


rex [oefes[os[os|oalar[os 

210.0} 210.0 | 290.0 | 250.0 | 210.0 | 170.0 | 120.0 
5x10° 0.069 | 0.069 | 0.102 | 0.094 | 0.086 | 0.077 | 0.063 | 0.050 

9 Moment | 234.0 | 290.0 | 380.0 | 330.0 | 330.0 | 290.0 | 200.0 | 140.0 
5x10° 0.108 | 0.098 | 0.148 | 0.126 | 0.139 | 0.137 | 0.107 | 0.090 


bx10° 0.155 | 0.243 | 0.214 | 0.166 | 0.177 | 0.180 | 0.155 | 0.119 
SSI fel fe 

ox10° 0.318 | 0.308 | 0.293 | 0.188 | 0.208 | 0.203 | 0.196 | 0.162 

6x10° 0.422 | 0.404 | 0.365 | 0.245 | 0.237 | 0.230 | 0.232 | 0.185 
eeeiceimmelcelelaaicals 

6x10? 0.546 | 0.489 | 0.513 | 0.285 | 0.277;:| 0.268 | 0.288 | 0.204 
7 ciceage 556.0 | 639.0 | 693.0 | 578.0 | 571.0 | 509.0 | 410.0 | 267.0 
x10" | 0.610 | 0.581 | 0.580 | 0.342 | 0.333 | 0.317 | 0.332 | 0.235 
Moment | 562.0 | 666.0 628.0 | 588.0 | 517.0 270.0 
6x10" 0.677 | 0.733 0.420 | 0.414 | 0.447 0.269 


Moment | 592.0 556.0 | 600.0 
x10° | 1.077 0.553 | 0.543 


6x10° | 1.136 0.668 
Moment | 597.0 

11 . 
x10> | 1.410 


Units: moment is inch-kips, curvature is 1/inch. 













TABLE 4.2 


Coordinates for linearized M-P-¢ diagrams for Type 2 column. 
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Axial load ratio, Pee 


6x10" 0.072 |0.069 |0.129 {0.079 {0.108 {0.099 |0.053 |0.057 
$x10° 0.118 {0.103 |0.220 |0.134 {0.151 | 0.151 {0.091 |0.086 
$x10° 0.227 |0.166 |0.268 |0.172 |0.183 {0.187 |0.132 |0.118 
6x10° 0.300 {0.196 |0.323 |0.292 {0.212 {0.228 |0.174 |0.144 
ox10° | 0.468 [0.285 |0.411 |0.356 |0.236 |0.262 | 0.219 |0.174 | 
$x10° 0.579 {0.363 | 0.516 |0.420 |0.268 |0.285 | 0.232 |0.190 
$x10> 0.731 |0.467 | 0.637 |0.470 |0.348 |0.300 |0.250 |0.213 
6x10° 0.862 |0.540 {0.707 {0.557 | 0.396 |0.356 | 0.285 | 0.260 
sh Siento ae Rl 
x10" | 0.947 | 0.630 |0.772 |0.670 |0.463 | 0.467 | 0.325 [0.304 
Sires <b secalr| 
x10" | 1.089 |0.691 Goole tO.ose_).0.006 
1] 


_ Point 
~s OY oO > Ww Po —_ 


©O 


Moment 500.0 296 .0 


3 


6x10 0.908 0.436 


Units: moment is inch-kips, curvature is 1/inch. 


TABLE 4.3 


Coordinates for linearized M-P-» diagrams for Type 3 column. 
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beyond this line the column and beam structures represented become 
unstable in a sidesway mode.The line mn will therefore be considered 
an instability failure condition. The line of instability corresponds 
to a point slightly above point 2 on the M-P-¢ curve. This type of 
instability occurs because equilibrium between external and internal 
forces is no longer possible due to the P-A moments and the reduction 
in member stiffness with increasing moments. 

In the balance of this report the type of failure described 
by Fig. 4.3 will be referred to as a "material failure" while that 
corresponding to Fig. 4.4 will be called a "stability failure”. 
Because the distinction between the two is difficult in some cases, 
the name stability failure has been reserved in this report for 
structures in which the moment at which the structure becomes un- 
stable is more than five percent below the moment capacity of the 
cross-section. 

The effect of slenderness and second-order effects is 
evident from the comparison of Figs. 4.3 and 4.4. For K=100 
the P-A moments account for about two percent of the failure 
moment for the short column and about 31 percent for the long column. 
For infinite beam restraint the values are 0.6 and 21 percent 


respectively. 


4.4 Normalized Subassemblage Charts 


One of the major difficulties in any attempt to develop a 
standard set of sway subassemblage charts for reinforced concrete 
columns is the wide variation in the M-P-¢ curves due to variations 
in geometrical properties, steel percentage, axial load and material 


properties. As a result, each particular cross-section has a unique 
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set of charts which cannot be applied to any other cross-section. 
This puts severe limitations on the use of the subassemblage concept 
for reinforced concrete structures since extensive calculations would 
have to be carried out for each case. 

An attempt was made to develop a normalized version of the 
M-P-¢ curves which would apply to a wide range of columns, but 


insufficient accuracy was attained for the methods tried. 


4.5 Deflections 

Few building codes specify limits on the deflection index 
A/1. In North America the National Building Code of Canada (32) 
is the only one to do so and the limit is set at 1/500 at service 
loads both for the story rotation and for the buildings as a whole. 
The same limit has been suggested by ACI Committee 435 "Allowable 


Deflections" >). 


When code prescribed values do not exist, values 
ranging from 1/300 to 1/1000 have been used in practice, depending 
on the type of building and the judgement of the engineer. 

A rigorous deflection analysis for concrete structures 
requires a considerable amount of computation and as a result approx- 
imate methods are usually preferred. 

However, the subassemblage analysis provides a better esti- 
mate if a realistic M-P-¢ curve is used. The deflections obtained 
from the subassemblage charts have been used for comparison with the 
deflection limits accepted in practice and to study the relationship 
between loading, deflections, geometry and mode of failure. 


Service load conditions were considered for two cases based 


on Equations (4.4) and (4.5) 
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U 


1.4D + 1.7L (4.4) 


U = 0.75 (1.4D + 1.7L + 1.7W) (4.5) 


where D,L and W are dead loads, live loads and wind loads, respect- 
ively. These equations are identical to Equations (9.1) and 
(9.2) in ACI 318-71. 

The two cases considered were: 

1. Maximum service load based on Equation (4.4) 

with L/D=5.0 
and 2. Minimum service load based on Equation (4.5) 
with L/D=0.2 

The service load was assumed to be given by Equation (4.6) 

Pe =D+L (4.6) 

From these relationships the axial service loads are 
equal to 0.61P and 0.92P for case 1 and 2, respectively. The 
corresponding load factors for lateral loads are 1.7 and 1.275. 
Thus the service lateral loads were taken as 1/1.7 and 1/1.275 times 
the maximum QL/M value for the ultimate load case. 

In figs. 4.5 through 4.14 the lateral load-deflection 
curves for ultimate loading conditions have been plotted for 
10 cases in addition to similar curves for the two axial loads 
corresponding to the maximum and minimum service load conditions. 
The service lateral loads have been obtained by dividing the lateral 
load by the appropriate load factors and the service load deflections 
were obtained from the appropriate load-deflection curve. The results 
are presented numerically in Table 4.4. The curves plotted in Figs. 
4.5 through 4.14 were chosen to study the effect that the shape of the 


cross-section, the slenderness and the beam restraint have on deflect- 
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ions. 

It can be seen from these figures that the shape of the cross- 
section has little effect on the deflections. The slenderness ratio 
is important for service load deflections and is also significant for 
ultimate deflections for strong beam-weak column combinations. As 
expected the deflections increase with reductions of the beam stiff- 
ness. This effect is more pronounced for material failures than for 
stability failures. 

All the deflection indices obtained in the analysis are 
shown graphically in Figs. 4.15 through 4.22. A total of 690 cases 
were studied, 217 of which were stability failures. A stability 
failure was assumed to occur when the moment at which the structure 
became unstable was more than five percent below the moment capacity 
of the column. 

Deflection indices at ultimate (Fig. 4.15, 4.18) showed 
wide variation with results ranging from about 1/700 to 1/30. 

The results for the two service load conditions 
(Fig. 4.16, 4.17) gave a somewhat more narrow distribution with 
values between 1/760 and 1/50 for case 2 and 1/820 to 1/70 for 
case 1. The service load deflections may vary from about 10 to 80 
percent of the ultimate deflections depending on the shape of the 
load deflection diagram and the service load condition. If the 
curve is very steep the service load deflections will reach the 
upper limit; in the case of large ultimate deflections the value 
will be closer to the lower limit. Material failures will normally 
have relatively small deflections while stability failures exhibit 


much larger deformations at ultimate. Typical examples are shown 
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in Figs. 4.5 through 4.14. 

Structures failing as a result of material failure were 
found to have service load deflections of about 1/3 to 2/3 of the 
ultimate deflection. If the loading condition corresponds to mini- 
mum service load the values will be in the lower range, say 1/3 
to 1/2, with weak beam-strong column combinations giving the higher 
values. For maximum service loads the values range from 1/2 to 2/3 
of the ultimate deflection. 

In Figs. 4.20 through 4.22 the computed values have See 
compared to allowable values of A/L of 1/500 and 1/250 for service 
load and ultimate load, respectively. For all three cases only a 
small percentage of the results falls below the specified limits. 
Thus, for the columns considered, deflections would likely be the 
governing design criterion rather than strength. 

Stability failures usually occur at relatively large 
deflections. Fig. 4.19 shows that no stability failures occurred 
for A/L less than 1/200. For allowable values below this figure 
the design would be governed by deflections or strength. This 
limit is used to set limits on the applicability of the PA type 
of analysis in Section 6.43. 

To investigate a possible relationship between deflection, 
loading and slenderness ratio the deflection index at ultimate was 
plotted against slenderness ratio for various values of axial load, 
Fig. 4.23 through 4.25. These plots show that the relationship bet- 
ween slenderness ratio and deflection index is close to being a 
straight line for relatively small values of kl/r. For higher 


values of k1/r the results become more scattered and the deflection 
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index reaches values which would normally be unacceptable in design. 
A straight line relationship between kl/r and (4/1), was therefore 
chosen such that it would give a reasonably good estimate of small 
deflection indices. 

The following empirical equation was taken to represent 
the deflection index at ultimate as a function of the slenderness 


ratio of the column and the axial load: 


ESV 2 Sees eerie Op. 4.7) 
D) (4.7) 


This equation has been plotted in Fig. 4.23 through 4.25 
for eae equal to 9.1, 0.4 and 0.7. 


4.6 Column Stiffness Analysis 


It was previously pointed out that the member stiffness 
is probably the most important parameter in any second-order analysis. 
A linear-elastic analysis will yield good results if the assumed stiff- 
ness is close to the effective elastic stiffness for the particular 
point in the load-deflection diagram being considered. 

Before determining an equivalent linear-elastic value of 
EI a desicision must be made as to what point on the load-deflection 
diagram should be used since the stiffness depends on the magnitude 
of the load and the loading history of the structure. 

If the coordinates of a point on the load-deflection curve 
are known, the equivalent stiffness may be computed from the following 
equations, derived from Fig. 4.2 using standard elastic formulae. 


The equivalent column stiffness is given by 
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me 
(EI). = ae, a (4.8) 
3(A - 3K, ) 
(EI), 
where Ky = aes (4.9) 
From Equation (3.25) 
3(EI), 
M. = KOM, = i 8 (4.10) 
KM 
Ke = | (4.11) 


The moment M is given by Equation (4.2). 

Fig. 4.26 shows a typical load-deflection diagram. If 
this structure had been loaded up to point B the elastic load-deforma- 
tion relationship corresponding to failure would be given by the line 
OB. If, on the other hand, the allowable or ultimate defTections 
limit the ultimate load on the structure to that corresponding to 
point A the ultimate deflection would be Ay rather than Ao given by the 
EI based on Tine OB. This overestimate will result in a corresponding 
overestimate of the P-A moments, in this idealized case of about 60 
percent. Similarly, if a point below A is chosen for computing the 
equivalent stiffness, an unsafe estimate of deflections and second-order 
moments will result. 

The discrepancies pointed out above will normally be severe 
only for structures exhibiting a very ductile behaviour. For a load- 
deformation characteristic such as that shown in Fig. 4.6 only minor 
errors would occur. 

It will be appreciated from the above discussion that choos- 


ing the correct point on the load-deflection diagram to define EI will 
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MATERIAL AND STABILITY FAILURES 





82 


SAYMNTIVA ALITIEVILS 
Sd9V4 UNOS NI SUVE *NWN109 GAIL YO4 SANTWA SSANSSILS 3AILN9449 §©OS*b “OIA 


zUt-qL “TW 
OL X 00*o8e 00 'Oh2 00 ‘002 00°09} 00 ‘021 00°08 00°Oh 
Vv 
¥ w 
v v 
eee: 
wwvy®y 3s ¥ 
¥wdy : 
vv w ¥ f 
ee 
v 
¥ 8 % 
vy Yv v ~- 


(8"O[) “ubz “Toy 


Z£°OL)*ubz * TOV 


00°O 


00°0 


00°OS1 00°00} 00°aS 


00 °00e 


A/9 





83 


OL X 00'082 


SSUNTIVA ALITIEVLS ONY TWIYALWW 


NWN109 TWuIdS YOs SANIWA SSIN44ILS 3AILO9II3 
| sul-ql *13 
00 'on2 00 ‘002 00°091 00 *021 00 ‘08 


S (SOL) uby “TOV 


LE fe_s9I4 


00 *Oh 


00°90 


00 °OQOT 00 °OS o0°O 


00°OSt 


00 *002 


A/O 





84 


Ol x 00'082 


SAYNTIVA ALT TELS 
NWN109 TWHIdS YO4 SANIWA SSANSAILS 3AIL094453 «= 2E*b “OIA 


ZuL-qL ‘I3 
00 ‘ah2 00 ‘002 00 "091 00 *021 00*o8 


@ 
fa) 
ge 
ar? 
y @ Qe e @ 
eae Ee 
2 @ & o @ 
ap 
- o 8 @ 
oe @ 
@ @ 
@ 
@ 


(8°OL) Ubz “Toy 


(Z°OL)*ubz *T9y 


00 °Oh 


00*O 


oo*Ost o00°Cot 00°OS 00°0 


00 °00e 


A/% 





' 


A en ey ee a ee ee 
‘OF x 06,085 00. 05 09 .0Os 0G G3! 00.08! 80 .68 


“ntedt .33 


WUJOD JABIG2 AOI 23UJAV 22IMAAIT2 BVETOSIW se.e OR 
23RUSEAR YTLSIAAT2 


66 


(c) 


Q— 





° q n 


FIG. 4.33 CONSTRUCTION OF LOAD-DEFLECTION CURVE 
FOR SINGLE STORY INTERIOR SUBASSEMBLAGE 
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depend on several factors and becomes a complex matter. The point corr- 
esponding to the maximum value of QL/M was therefore chosen to be used 
in the calculations. This will always give conservative estimates of 
both first and second-order deflections. 

The computed effective linear-elastic stiffness values have 
been plotted against slenderness ratio in the form of scatter diagrams 
in Figs. 4.27 through 4.32. A comparison with values computed from ACI 
318-71 equations (10-7) and (10-8) shows that only a small number of the 
computed results are below the ACI values, and that for a large number of 
cases the code equations give extremely conservative results. The few 
cases where the effective EI is smaller than that computed from the ACI 
equations corresponded to low values of PP: For erie < 0.2 the most 
conservative of the ACI equations gives the best estimate of EI. 

In the case of stability failures the columns tend to have ra- 
ther high stiffness because the curvature at failure is relatively small 
and thus the effect of cracking is less severe. High EI values were also 
found to occur for material failures when the axial load was fairly close 
to the balanced load. 

It should be noted that the ACI Moment Magnifier procedure treats 
stability failures as material failures. To do this it is necessary for 
the EI value for use in the ACI procedure to be based on the point of in- 
tersection of the load deflection curves with the sloping curve represent- 
ing the cross-sectional strength. Thus, in the case shown in Fig. 4.4, 
the effective EI corresponding to the maximum value of QL/M (points a- 
long line mn) will be considerably larger than the EI corresponding to 


the eventual failure of the section (points along line ab). 


4.7 Use of Sway Subassemblage Diagrams in Frame Analysis 
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The response of a story subjected to lateral loads and gravi- 
ty loads can be determined by superposition of the response of each in- 
dividual column in the story. In order to establish the column resp- 
onse the story is subdivided into a number of subassemblages as out- 
lined in Section 2.6 and the load-deformation characteristics of each 
subassemblage is determined using the restrained column curves. 

A typical subassemblage of an interior column in a single 
story frame is shown in Fig. 4.33(a). The restraining action provided 
by the beams is K 


6M for beam AB and K OM, for beam BC. The restrained 


1 2 
column curves relevant to this subassemblage are shown in Fig. 4.33(b). 
Line o-d | is the restrained column curve for a column restrained by 

both beams, line o-f : and o-e | are for the same column restrained 

only by beam AB or only be beam BC, respectively. The lines m-n and 

p-q correspond to the formation of plastic hinges in beam BC and beam 

AB respectively. 

The load-deflection curve for the subassemblage has been con- 
structed in Fig. 4.33(c). taierhike the load-deformation response is 
given by the line Oa which is parallel to 0a Wwrid. 4.5000). At. a’ 

a plastic hinge forms in beam AB and a further increase in load is re- 
sisted only by the beam BC. The line ab is parallel to the line be 
in Fig. 4.33(b). At b a plastic hinge is formed in the beam BC as 
well and the structure becomes unstable. Thus the ultimate load capa- 
city of this subassemblage is given by point b. 

If either one or both of the beams had a moment capacity equal 
to or greater than that of the column the load-deflection curve would 
reach the upper inclined line. If neither beam developed hinges the load 


deflection curve would be given by oa d ,» and if only the weaker beam 
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developed a hinge the line ab would be replaced by bbe: Gi 

When the load-deflection curve has been constructed for all 
subassemblages in a story the total shear resistance may be determined 
by a graphical super-position of the individual load-deflection 
curves. 

In a multi-story frame the overturning moment produced by the 
column above must be considered (see Section 3.5.1 and Fig. 3.11). 
If restrained column curves based on the model in Fig. 4.2 are to be 
used for a multi-story frame a modification must be made to the beam 
stiffness K. From Equations (3.23), (3.25) and (4.1) the following 
two expressions for the joint rotation 6 are obtained for a single- 


story and a multi-story frame respectively: 





M, 
— (4.12) 
u 
2M M 
e=_t=— 2 (4.13) 
KM (ym 
a ou 


Thus, for a multi-story structure the sway subassemblage load-deflect- 
ion curves should be based on restrained column curves corresponding 
to a beam stiffness of one half of the actual beam stiffness. 
References 17,34 and 35 present design examples and design 
aids for the design of multi-story steel structures by the subass- 
emblage method. These design curves differ from those presented here 
in that they are based on the assumption that an upper column frames 
into the joint and thus, if the curves are used for a single story 
frame the curve corresponding to a beam twice as stiff as the actual 


beam should be used to determine the load-deflection response. 
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CHAPTER V 
EFFECTIVE STIFFNESS OF REINFORCED CONCRETE BEAMS 


5.1 Introduction 

No firmly; established method exists for computing the effective 
stiffness of reinforced concrete beams. Various approximate methods 
are used in practice, some of which sacrifice accuracy for simplicity. 

However, a reasonable accuracy in the estimate of the restrain- 
ing action of the beams is essential in any second-order analysis and 
some guidelines are therefore required for its computation. 

In this chapter the effective stiffness of reinforced concrete 
beams will be studied for various loading conditions and loading in- 
tensities. An attempt will be made to establish some simple rules 


which may be used in second-order analysis. 


5.2 Method of Analysis 


A cracked reinforced concrete beam behaves as a beam having a 
varying moment of inertia along its length. The effective moment of 
inertia at any section of the beam is a function of the bending moment 
at the section and may also depend on the loading history of the struct- 
ure. 

To account for the effect of the change in bending moment along the 
beam the beam was divided into a number of small segments such that 
the bending moment in each segment is approximately constant. Thus it 
is assumed that the distribution of cracks is uniform throughout each 
segment. 


If the moment-curvature relationship of the cross-section is 
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known the stiffness of a segment having a certain bending moment is 
easily obtained as the slope of the moment-curvature diagram at the 
point being considered. To simplify the computation of the slope, the 
M-o¢ Curve was approximated by a number of straight lines as prev- 
iously outlined in Section 3.4.4. 

The moment capacity of the cross-section was assumed to 
be reached when the tension steel yielded, thus ignoring the effect 
of strain hardening. This assumption gave failure moments in close 
agreement with values obtained from standard design formulae for 
computing ultimate moments. 

The moment in each segment was assumed to be the mean of 
the moments at the ends of the segment and from the M-o curve the 
stiffness, EI, corresponding to the mean moment was determined. The 
flexibility coefficents were obtained by loading the beam with the 
M/EI diagram corresponding to unit end moments and applying the con- 
jugate beam theorems to compute the resulting end rotations. The 
rotational stiffness coefficients were then established by inverting 
the flexibility matrix. 

It should be noted that the distribution of moments was a 
function of the distribution of stiffnesses which in turn were a 
function of the distribution of moments. An iteration procedure was © 
required to reach a distribution of moments and stiffnesses which 
were compatible with the applied loads. The rotational stiffnesses 
were calculated based on the final distribution of stiffnesses. 

It was found that the redistribution of moments did not 
cause significant changes in the stiffness values. This could also 


be anticipated by considering the idealized shape of the M-¢ curve 
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since only elements which initially had a moment value close to 
points of discontinuity would have its stiffness affected by a 
change in moment. The idealized M-¢ curves have only one or two 
points of discontinuity and as a result only the stiffnesses of a 
few elements are affected by the redistribution of moments. 
Effective stiffness values were computed for three 
loading conditions: 
1. gravity load only. 
2. lateral load only. 
3. gravity load and lateral load combined. 
Each case was investigated for various values of load intensity. 
The load intensity may conveniently be expressed as a 
function of the end segment moment capacity. In the case of gravity 
load the beam was assumed to have fixed ends and the nominal 
fixed-end moment Mom may be expressed as 
2 


wl) 
em = 2 = 7M, (5.1) 


where w is the uniformly distributed load on the beam, us is the 


M 


length of the beam, My is the moment capacity of the end segment 

and n is a coefficient which may have values between zero and one. 
The end moments, Mow? caused by lateral loads can be 

written as 

sm = My (5.2) 

where » is a coefficient varying between zero and one. 


M 


For combined loads the applied end moment at the most 
heavily loaded end is given by 


(Ma, + M =- (n +u M,, (5.3) 


GM sm) 
and at the other end 
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(Moy = Mow) (n 7 uM (5.4) 
assuming counterclockwise moments to be positive and the lateral 


load to be applied from left to right. 


5.3 Effective Stiffness of T-Beams and Rectangular Slabs 


Two different shapes of cross-section were considered in 
the analysis, a T-section and a rectangular flat slab section 
(Figs. 5.1 and 5.2). 

The T-section had a total depth of 22.5 inches and a web 
width of 13 inches. The flange was 4.5 inches deep and 78 inches 
wide. The effective span was 25 feet. The slab was 12 inches wide, 
7.5 inches deep and had an effective span of 18 feet. 

Both members were symmetrically reinforced about the 
centre line and were divided into two end sections and a middle 
section. All three sections had equal amounts of top and bottom 
reinforcement, the ends being more heavily reinforced than the 
middie. The details of the cross-sections are shown in Fig. 5.1 
and Fig. 5.2. 

Typical M-¢ diagrams obtained for the T-section and the 
Slab are shown in Figs. 5.3 through 5.5. Fig. 5.3 shows the response 
of the T-section when loaded such that the flange is in compression 
while in Fig. 5.4 the flange is loaded in tension. The latter case 
shows a large and sudden drop in moment capacity as cracks develop in 
the flange. However, this is an artificial equilibrium condition 
which shows up in the computations because the moment capacity is 
computed from predetermined increments in the curvature. The transi- 


tion line from the point of cracking to the point TY was used in the 
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SECTION A-A AND C-C 
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FIG. 5.2 SLAB FOR STIFFNESS ANALYSIS 
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computations. Fig. 5.3 shows a similar decrease in stiffness after 
the initiation of cracking, but the transfer is much smoother. The 
M-¢ response of the slab shown in Fig. 5.5 is similar to that of 
the T-section with the flange in compression. 

A few trial calculations were carried out to investigate 
what length of segments should be used to attain sufficient accuracy. 
It was found that for a given length of segment the accuracy was 
dependent upon the type of loading, and the largest inaccuracies 
seemed to occur for combined gravity and lateral loads. If high 
accuracy is sought a fairly short segment length should be used. 

In these calculations 40 segments were used for the 
T-beam and 72 segments for the slab. 

The results of the calculations have been plotted in 
Fig. 5.6 for the T beam section and Fig. 5.7 for the slab section. 
Considering the members under pure gravity load (Figs. 5.6a and 
5.7a) shows clearly the reduction in stiffness as the load intensity 
is increased and the cracked zones extend. 

The T-beam is seen to keep 95 percent of its uncracked, untransformed 
stiffness E. IG? for values of n below 0.5. For higher loads the 
stiffness is reduced rapidly and for n = 0.9 it has dropped to 46 
percent of the uncracked stiffness. 

The effects of cracking become evident at fairly low load 
levels for the slab and once the first cracks have opened up there 
is a rapid decrease in stiffness. However, as the load intensity is 
increased the rate of stiffness reduction is reduced. When the moment 
at the ends has reached 90 percent of the capacity of the end segment 


the stiffness has dropped to 53 percent of the uncracked stiffness. 
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In a second-order analysis the beam will be loaded 
according to load case 3. The variation in stiffness for various 
combinations of fixed end moment and sway moment is shown in Fig. 
5.6(c),(d) and Fig. 5.7(c),(d) for two values of load intensity. 
The plots show only small variations in stiffness with variations 
in the ratio n/u. For very small values of n/u, say below 1.0, 
corresponding to high lateral load moments and small gravity load 
moments it would be more appropriate to use value from load case 2, 
i.e. sway moments only, as shown in Fig. 5.6(b) and Fig. 5.7(b). 

As one would expect, reducing the load intensity (ntu) results in 
higher stiffness values, although for the slab this increase is 
negligible. In practice one would expect to find values of (ntu) 
in the range 0.75 to 0.90 and the results presented herein suggest 
that values of (EI), between 0.5 and 0.6 times El, may be used for 
this type of loading condition if n/u is greater than 1.0. 

The steel percentage will of course affect the stiffness 
of the member to some extent. This maybe taken into account by 
expressing (EI), as 


(EQ) em beled. tly) (5.5) 


g 
where a is a coefficient, n is the modular ratio of steel and con- 
crete and iz is the moment of inertia of the reinforcement about the 
neutral axis of the gross concrete section. 

For load case 3 and n/u greater than 1.0 the value of the 
coefficient ‘a may be taken as 0.2 for the T-beam and 0.15 for the 
slab. For a modular ratio of 8 this gives the stiffness values 


0.496E 1 for the T-beam and 0.53E 1, for the slab. 


g 
The effective moment of inertia of T-sections is often 
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expressed in terms of the moment of inertia of the web alone, I 
In this case the ratio L/ly is equal to 2.0 and the effective (EI), 
is therefore between 1.0 and 1.2 times El: This is somewhat 


36) 


lower than a commonly accepted ( assumption that (EI), =2(E 1) 


and higher than a method assuming the effective flange width equal 


to twice the web width (37) 


2 2(E.1,) is equivalent to an effective 
flange width of about three times the web width. 

The stiffness parameter K, used in the sway subassemblage 
analysis, may be computed for these beams by using Equation (4.11). 
The moment capacity of the column, Mi will be assumed to vary between 
828 k-in and 300 k-in (Type 1 column, see Fig. 4.1 and Table 4.1). 
The T-section will then have K values varying from about 500 to 
about 2000 before cracking and only about one half of these values 
after cracking. For the slab the values vary between 50 and 


140 before cracking and again about one half of these values after 


cracking. 
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CHAPTER VI 
ANALYSIS OF MULTI-STORY FRAMES FOR STABILITY EFFECTS 


6.1 Introduction 


There are two types of stability failures that may be 
encountered in a frame analysis. One is known as "bifurcation of 


equilibrium," or buckling, and occurs when the applied axial 

load reaches the critical buckling load. The other type is referred 
to as “instability through disturbance of equilibrium" and occurs 
because equilibrium between external and internal forces cannot be 
achieved due to such things as imperfections and reduction in stiff- 
ness. This type of instability occurs for an axial load smaller 

than the bifurcation load of the member. In sway frames where the 
members exhibit a load-deformation response as shown in Fig. 6.13(a), 
sidesway instability will occur before the ultimate moment has been 
reached. When instability occurs, a small increment in the lateral 
load produces additional PA moments, which in turn cause additional 
deflections. Because the applied moment is a function of the lateral 
deflection and because the stiffness decreases with increasing 
moments, the required internal moment resistance is not achieved. 

The PA effect is, therefore, a key parameter which must be considered 
in the analysis of sway permitted frames. 

Modern building codes attempt to predict stability failures 
by means of simple approximate methods. The accuracy of these methods 
may be quite good in some cases, while in others they are highly in- 
accurate but generally conservative. 


A rigorous stability analysis of reinforced concrete frames 
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is a rather complicated matter due to such things as the non-linear 
load-deformation relationships of concrete columns and beams and the 
effect of the steel percentage and axial loads on the member stiffness. 
A general method which would give good accuracy in all cases would 
have to consider the effect of these variables, resulting in a very 
complicated procedure. 

The subassemblage procedure provides a good method for 
predicting instability and the results from such an analysis will be 
used as a basis for comparison when investigating the applicability of 


current procedures for second-order analysis. 


6.2 Sway Preventing Action in Frames 


Prior to embarking on the design of columns in a frame 
by traditional means it is necessary to determine whether the frame 
is braced or unbraced since the behaviour of these two types of frames 
differs greatly. The problem is how much lateral restraint is re- 
quired to allow a frame to be designed as a braced frame. 

The ACI Code Commentary <=) states that the bracing elements 
(shear walls, trusses, etc.) should have a lateral stiffness of at least 
six times the lateral stiffnesses of all the columns in the story being 
considered. This requirement may not apply equally well to all struct- 
ures, however. 

A more rational evaluation of the sway preventing action in 


(38,39) The analytical model used by 


frames has been presented by Lay 
Lay is shown in Fig. 6.1. The end conditions of the column are repres- 
ented by the rotational stiffness coefficients Kn and kes and the trans- 


lational stiffness coefficient Ko. The restraining actions at the end 
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FIG. 6.1 LAY'S MODEL FOR INVESTIGATING SWAY PREVENTING ACTION IN FRAMES 
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of the member can be written as 


™m 
— 


My = ky Ts (6.1) 

My = ky =P 6p (6.2) 

Qn = k, 7 (*) (6.3) 
The equilibrium equation is 

PA = MatMa+QpL (6.4) 


My and Mp may be determined from the slope-deflection equations 
and substituted into Equation 6.4. 


The solution is then: 


2+s(1-c) (Ls oT 
Ses D k.* k 
k. = (5—) - s(l+c) A “B 
s Pel ar $3 (1-c2) (6.5) 
ky kg Kgke 


where gar is the buckling load, s and c are the standard slope-deflect- 
ion coefficients modified for axial load and are given by Equations 


(6.6) and (6.7) 





ein Ton OEM (4) (6.6) 
tans->5 
eho hs 
c=ucscu - 1 (6.7) 
1 - cot u 
u = TY P/PE . (6.8) 
2 
_ a7 EI 
Lay (39) suggested three ways to use Equation (6.5) to investi- 


gate the sway stiffness necessary to ensure sway prevented action in 
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frames: 
1. An exact evaluation of k. using Equation (6.5) 
2. If the numerator of Equation (6.5) is positive and 
Kn and kp are both positive, which is always the case 
if P/PE < 1.0, the lateral stiffness necessary to allow 
the column to reach its no-sway buckling load may be 


approximated by: 


ONT A Oe (6.10) 
cr 


This will give a conservative value of the critical kK. provided the 
rotational end restraints are positive. Negative end restraints will 
rarely occur in ordinary building frames due to the low axial load 
usually carried by the beams. Thus, this will be conservative provided 


that the numerator is positive, i.e. 


2+ s(I-c)( + aD = '0 (6.11) 


> 


If the required lateral restraint exists according to Equa- 
tion (6.5) or Equations (6.10) and (6.11) the effective length factor may 
be determined from the nomograph for braced frames. When doing so, 
Trahair's method (40) (Equation (6.21), Section 6.3.2) should be used to 
determine the parameter y. Because the frame is braced the effective 
lengths will be less than the actual lengths. 

3. The buckling load may be approximated by Equation (6.10) 
without any further checks if the effective length is not 
taken less than the actual length. 

The derivation of these stiffness criteria was based on a 


linear-elastic load-deformation relationship. 
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Telwar and Cohn (41) have suggested four bracing criteria 
for tall buildings containing shear walls: 

1. To limit the deflections at service load to a value 
of 1/500th of the story height they proposed that: 


p 
700 
re ON (6.12) 


where - is the critical load for the entire structure and Q and P 
are lateral and vertical service loads. 
2. A stability criterion to ensure that the frame could 


be designed for gravity loads only was the basis of Equation (6.13): 


P 

+= 2 90 (6.13) 
3. They also proposed that 

ie 


Por 


t= 
lv 


0.5 (6.14) 


where Bo is the critical load for the wall standing on its own. 
4. The fourth criterion was intended to prevent excessive 
moments due to unsymmetrical loading and for a pseudoelastic frame 


was expressed by Equation (6.15) 


12k 0 f 


Pop = Max| (7) (==) ae (6.15) 


c 

where i is the story under consideration (counting from the top), 
h is the story height, k. is the EI/h for the column in the story 
carrying the greatest shear and af and O° are the shears on the 
frame and the column, respectively. The third and fourth criteria 


are based on assumptions about desirable behaviour and may fail to 
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recognize the entire spectrum from fully braced structures through 
structures free to sway. 

One or more of the criteria presented in this section may 
be used to evaluate the sway preventing action of the bracing 
elements in a frame before a method of analysis is decided upon. 

A further study of the sway preventing action in frames 


is presented in Section 6.4.2. 


6.3 Effective Length Factor Method 
6.3.1 Introduction 


Traditionally the effect of frame action has been 
accounted for in column design by means of effective length factors. 
This method has been described in Section 2.2. In this section the 
results obtained from this method will be compared to those obtained 
in the subassemblage analysis. Three factors are of prime interest; 
the effective length factor itself, the magnitude of the amplified 
moment and the mode of failure. 

6.3.2 Some Remarks on the Effective Length Concept 


In both the ACI column design method‘! and the CEB column 
(4) 


design method there are two very important parameters involved in 
determining the second-order effects, namely the stiffness, EI (or 
curvature in CEB), and the effective length factor k. The latter is 

a function of EI since it depends on the ratio of the sum of the 
stiffnesses of the columns framing into a joint to the sum of stiff- 


nesses of all the beams framing into the same joint. This may be 


expressed as 
Z(EI/1). 


p = TEI/T), (6.16) 
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This expression is an approximation to enable the effective length 
nomographs to be used for irregular frame lay-outs. 


(39) has investigated the validity of this approximation 


Lay 
and from an elastic slope-deflection analysis he developed the follow- 


ing expression for the relative stiffness parameter 


(EI/1) 4. ie 
oo : 
Bq) + 2/59 Fp), 


where (EI/1) 4. is the stiffness of the column being designed, 


du _ sQ-¢) 65. braced frames (6. 18a) 
Jo 2 
su SUC) Foy sway frames (6.18b) 


s and c are the slope-deflection stiffness coefficients including 
axial load effects as defined in Section 6.2. 

Once the end stiffness ratios are computed the effective 
length factor, k, can be obtained from nomographs presented in 
Reference 2 or from approximate equations such as those given by 
Furlong (72) or in the British Standard Code of Practice cp110'43) 

An examination of the nomographs shows that the value of k is 
reduced as the value of » is reduced, which in turn implies that the 
critical load of the column is increased. Thus it follows that the 
discrepancy between Equation (6.16) and Equation (6.17) will result 
in an unsafe estimate only when 

eae (6.19) 
which will occur when 
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where the subscript mc refers to the upper column framing into the 
joint. 

From Equation (6.18) it can be seen that this can only occur 
for braced frames and substitution of Equation (6.6.) into Equation 
(6.18a) shows that a second condition is that u (Equation 6.8) must 
be greater than 7. Thus, Equation (6.16) will lead to a safe estimate 
of k in braced frames if P < Pee 


To simplify the use of Equation (6.17) Lay suggests that an 


(41) 


approximate method developed by Trahair may be used. This method 


approximates Equation (6.17) by the following expression 


(FF) 
aes 1‘ column being designed 
Py El (6.21) 


2(1-5-) (7-) all other members 
E 


The term (1-P/P -) approximates the effect of axial load on the member 
stiffness and will normally be significant only for the columns since 
the beams usually have negligible axial load. 

Although Equation (6.16) will give conservative results 
for sway frames Lay (39) suggests that the modified Equation (6.21) 
be used since there are considerable economic advantages to be gained 
from it. 

Rosenblueth'”?, MacGregor, Breen and Pfrang'?Jand, Spring- 
field and Adams (45) have pointed out shortcomings in the effective 
length concept in dealing with sway frames where the columns in a parti- 
cular story have widely varying effective length factors. The extreme 
case of a pin-ended column supported by a sway frame occurs frequently. 
The nomographs suggest that the pin-ended column will not be able to 
carry any vertical load since k = ~ because sway can occur, and there- 


fore ae = 0 for the hinged member. Some modification of the method 


111 



















att niawloo 1u say of 218to7 on tafswedue of srg - 
“into 


Sens neg nso 3¢ (81.3). notésupd mov? yi 
pry ee 
o436u03 to mets hutttz¢ dvz bre zomext bsssyd. 70? 7. 


~~ 


rt ADD bnooee2 4 ~ aft owore (68 ay 


f (al iortoupa ~gudT .« men? Y9tss1D od 
.4t.> 9 W-zomen? beosid ot 4 te 
J ‘ } to seu ene ytriqmt2 oT See 
~ tel 
; een boqolfeval orl? yn a2 smi xo7%¢ 16 
: u ET af yd yi .B) 1O! ups 28d emt xo 
[3 
fons 
: ‘cs i ‘a teed) has ae | ‘ é 
= nathan nay oe i * we = ¥ F 
yon vite TTS (=e) e-f) 
~ = ‘ hy Cc é 


. 7, % te o ey 
bi xs atis sft extemal xoigge ig “V4. wits T. 


Ps 


“_. ’ id 3 ‘ee 
[fo Sag9rttheate sa A eel ftw brs zzo0 se 


ce 
Ns comet one 
: ee a 
wrap : i bs attonooe .sfdswsh tane9 a6 sions 92 aC 


- 


A ogy} s Need. TORO TESEM Nasoutdition ° 


. ; ease ae 
2p1 imoos xeric Juo! batatog sven ta 6), mada & tie bie 


At enmulo2 “SNS SIS ZomsyT Yewe dd iw sitet nt tq: phe 


ate vtxe siT .2y0to87 apne! evita tts: engib agile beahaees: 


anewpe3 e490. Smet? YEW? 5 xa besroqque maton haase 


112 


must therefore be considered where yw is very large for both ends of 
the member or where yw varies widely between columns in a frame. 

The 1971 ACI Code attempts to treat this case by replacing 
the P/P an term in the moment magnifier equation by the term zP/EP 
whose variations in effective length is recognized by treating the 


entire story as a unit (193), 


6.3.3 Comparison of Results from the Moment Magnifier Method and the 
Subassemblage Analysis 


The accepted practical procedure for determining the effect- 


(2,46,16,34)  ntioned 


ive length factor, k, is to use the nomographs 
in Section 6.3.2. The magnified moment, according to ACI 318-71, 
is then given by 
1- Ta se, (k1) 
c 
where Mo is the magnified moment, Mo is the first-order moment and 
Ch is an equivalent moment factor. 
Equation (6.22) may be rearranged to give an expression for 
the effective length factor k, and using the subassemblage model 
in Fig. 4.2 it can be written as: 
2 
jee. re ; d mesa 4 (6.23) 
Thus, the results from the subassemblage analysis may be 
used to compute the effective length factor which would have to be 
used in Equation (6.22) to give the correct second-order moment. 
The quantity y needed to find the value of k from the nomo- 
graphs is defined by Equation (6.16) with the denominator of the 


equation determined from Equation (4.10). The parameter (EI). was 
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calculated according to Section 4.5. 

In this case there is no upper column framing into the 
joint and thus Equation (6.16) is in agreement with Trahair's 
method, Equation (6.21). 

The quantity Ch will be taken as 1.0 in this derivation 
according to ACI 318-71 rather than the value of 0.85 as given by an 
elastic analysis for the sway case. This gives values of k somewhat 
below the correct value. Ch is a function of the first-order deflec- 
tion and since this is not known due to the non-linearity of the load- 
deflection curve, a precise value of C, is difficult to obtain for 
reinforced concrete columns. 

In Fig. 6.2 the effective length factors obtained from the 
nomographs and from the subassemblage analysis have been plotted against 
relative stiffness. The two curves show fairly good agreement, with 
the nomographs giving the higher values. The maximum percentage dis- 
crepancy between the two curves is about seven percent and occurs for 
low values of k. This percentage decreases with increasing k and for 
values of k about 3.5 the error is only one percent. If the correct 
value of C. had been used, the effect would have been to bring the 
two curves closer together, which indicates that the nomographs 
predict quite well the length of the equivalent pin-ended column. 

Figures 6.3 through 6.10 show the relationship between the 
short column interaction diagram (1/h=0), an interaction diagram 
for material or stability failure, taken as Point 1 in Fig. 6.13 (a) 
or (b) and a design interaction diagram based on the moment-magnifier 
procedure using Eqn.(6.22)and either the effective EI from Eqn(4.8) or 
the EI from ACI Eqn.(10-7). Results are plotted for values of end 
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6.3  LOAD-MOMENT RELATIONSHIP AT ULTIMATE 
OBTAINED FROM THE SUBASSEMBLAGE AND 
MOMENT MAGNIFIER METHODS OF ANALYSIS 
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restraint Ky = o and K a 10’ k-in, a value close to the lowest 
practical limit for beam stiffness, and various slenderness ratios. 
All plots are for a tied column with bars in two faces (Type 1 column). 
For 1/h greater than about 20 the failure line lies inside 
the interaction diagram indicating stability failures. Interaction 
diagrams based on the moment magnifier equation using the effective 
EI correspond reasonably well to the failure interaction diagram as 
they should,since the effective EI values were derived from Point 1 
in Fig. 6.13(a) or (b). If design were based on the use of this 
interaction diagram, the design would require less steel than 
actually necessary. If, on the other hand, design were based on the 
moment magnifier fi the ACI EI value, the design interaction diagram 
lies outside the cross-sectional interaction diagram indicating that 
more reinforcement is required than actually needed. Thus the ACI 
Code procedure predicts the second order moments well for short 
columns (Figs. 6.3 and 6.7) but tends to overestimate the moments for 


slender columns. 


6.4 Iterative P-A Procedure 
6.4.1 General Principles 

The basic principles of the P-A method were outlined in 
Section 2.4. This section will be devoted to a more detailed investi- 
gation of the method and its application in the design of building 
frames (47) 

Fig. 6.11 shows a diagram of a column in the i th 

story of a frame carrying both lateral and vertical loads. The deflec- 


tion A is computed ignoring the effect of Pat From equilibrium 


requirements the additional shears produced by the axial load acting 
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FIG. 6.11 SWAY FORCES DUE TO VERTICAL LOADS 
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through a deflection A may be written as 
' zP 
n- T '4ne1 - 4) Mh-#4) 
where Vv, = additional shear in story n due to axial loads, 
zP = sum of column axial load in story n, . = height 
of story n, 


A ‘7 deflection of floor level n+ 1 and n respectively. 


ntl, 


The net additional sway force at level n due to axial loads is given by 


Qn = V - V (6.25) 


This force is now added to the applied lateral load at level n and a 
first order analysis is again carried out. The deflections obtained 
in that analysis are substituted into Eqn. (6.24) to give a new set 
of story shears which in turn are used to compute a new set of addi- 
tional sway forces. The process is repeated until the moments or 
deflections in one cylce are only slightly larger than those of the 
previous cycle. 

The rate of convergence of the iteration process could 
be used as an indication as to the possibility of a stability 
failure occurring. Some results from the P-A analysis are shown in 
Tables 6.1 and 6.2, where A. is the deflection in the i th cycle 
and Ae is the final second-order deflection as obtained in the sub- 
assemblage analysis. The short column interaction diagrams and the 
failure curves for these columns are shown in Figs. 6.3 through 6.10. 
The iteration process is seen to converge slowly for stability failures 
and it appears that if the process converges within five iterations 


there is no danger of a stability failure occurring. 
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TABLE 6.1 


Results from the PA iteration procedure. 


ACI Eqn. (10-7) 
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Tables 6.1 and 6.2 also show how the deflections and the 
rate of convergence is affected by using EI values computed from ACI 
Equation (10-7). Because the Code equation underestimates the stiff- 
ness in most cases the deflections are overestimated and the ratio 
Ac/A; is therefore less than one. The rate of convergence is seen 
to decrease with increasing Pm or and as a result, using a smaller 
value of EI will reduce the rate of convergence. 

6.4.2 Direct Solution 

The process described in Section 6.4.1 may be expressed in 
a more convenient form. Let Qu and Ay represent the applied lateral 
load and the corresponding first order deflection, respectively. 
Also, let the appropriate axial load be P, the deflection caused 
by a unit lateral load be K. and let Q. (i=2,3,--~-) be the sum of 
the applied and additional lateral loads in the i th cycle. 

Then the iteration process may be expressed in the fol low- 


ing manner: 


lst iteration: A,=k,Q, (6.26) 
ALE, nd Sihtes p 
2nd iteration: do=k Qo=k.Q, (147k, ) (6.27) 
3rd iteration: A,=k-Q.=k_Q, (14k +(2)?k2) (6.28) 
“acs S aueS tle awl Soe LS i 
and the general term for the i th iteration is 
i-2 . i-l. 
. PoosaPaece eto. Pe i-l 
Ai=k.Q) 1+7k.+ (7) ketenes) k. +(7) k. (6.29) 


This is a geometric series which converges if + k.<1.0, and in that 


case the sum of the infinite series is 





: a 1 
1-~k 1 
1l°s ler mg 


128 




























Ss 
BSI ie kee 
eft bas znotiosfieb omg wor ware gels $ia tin { 4 esfdeT ” a 
[DA mov} baduqinoo zoutsy 13 paren yd vayatts et soneprevnes to ster 
-Ytife efi estantyeswsbAy avid supe obod sat gausosd - (1-01) note sup3 
“otis st brs bet emtsearsve a5 amor $9 Piast 2o262 d20m nt 228n- 
. sbvanaeel to stsx fT ono dem? e2et groterend at ing , 
selfeme s onfau .dfueex 6 26 Ons, 9 pntamanoat Agtw sesetosh oF - 
Sonepyevie2 to stay ohh Bouboy Et fy Id t0. oul sv 


notsuloe soext0 $0.8 


a) 
at beersvaxe od yom [.4.38 nottoee vt bsdivgesb senoe'n adT : 


fereate! betfieas of) $a9a97i97 1? bith p ta »ryor dnetasvags enon § 


aaa 
(fevttseqee .notsoal tsb ywbve raw pat bnogee709 ony bn beof 7 


be2us3: hot too Mabsond 74 st anol [arks adel qoregs ent set veal 


10 muz atit-ad (=, EST) 0 af bis r, a} baat ‘fenstef Show 6 ec 
stove dt t -onbint aback Teveter Fehots bobs bre “bat faqs ad 


wot fo* en? at bszeengx9 od yen 2899079 notisredt eds not. “as ae | 


~~ 


oe 


—. 
fe’D 


(BS .3) ; Paap suites gal 
(0.a) (APT) Ota Peet. inatenatt bas 9) = 
——, Soe a eis J aye ’ 
Cr+ te “£}.,0 al Naph seiscaamen ves oe 
at nofdexatt @ t - ro? erat. isranoe oft 6 


(28.2) [t- . Fé + pe , a " dail ‘a 
nf o's 


jdt wt bes Ot +? ents ie in ooeaa 


iat ~ % pSy5 
(o3.0) ‘3M 












0.5 
Failure M =4660 
0.4 ars for =0.4 


Pa Raed yi 
Exact EI 
~*) 


0.3 Eqn. 6.31 


200 400 600 800 1000 1200 
MOMENT K-in 


FIG. 6.12 P-A ANALYSIS COMPARED TO SUBASSEMBLAGE ANALYSIS 


SG" £3 
K. = k-in. 


es! 









Oddh= M 
gq 
AO: 10} 


of x 
bofltem A*9 AS 


3 + 
ale g: S25 -13 






“oost «600or «= 008 S(t 


212¥ JAMA yon savazengue OT GamiASMaD: penne ven ) 
= vein 23 - r 


ann bf Gate 9s “a5 











P, 





oe f ans 
' i : ‘2 4 * al _— _ q 
a hg a “ie » ee: . 7 . ? 
eet Sg thn i Se as 
es oe ty Stree 
Pi a ce 


130 


And since kK.Q3=4) the process converges to the final second-order 


deflection Ay. 


Ao = rere, (6.31) 


Equation (6.31) is identical to the equations proposed by 


( (10) (7) 


Fey 9) Parme and Goldberg and shows that the second-order 
deflections may be computed directly from the results of a first- 
order analysis in the case of a straight line moment-curvature 
relationship. If a non-linear M-P-¢ curve is approximated by a number 
of straight lines, Equation (6.31) may still be used if the lateral 
load is applied in increments to allow for appropriate adjustments 
of the EI value with increasing moment. 

The obvious advantage of using Equation (6.31) rather than 
the actual iteration process is that only two first order analyses 
are required to obtain the second-order moments and forces in an 
elastic structure. 

The accuracy of Equation (6.31) may be studied by considering 
Fig. 6.12. This figure shows interaction diagrams obtained from: 
a) the actual iteration process, b) using Equation (6.31) and c) the 
subassemblage analysis. The iteration process was carried out using 
both the effective EI and ACI equation (10-7). The results from 
Equation (6.31) were obtained using the effective EI. 

It can be seen that the use of Equation (6.31), with the 
correct EI, gives results identical to those obtained in the sub- 


assemblage analysis while the iteration procedure gives smaller 


moments. The moments from the P-A analysis were lower because the 
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P-A calculation converges on the correct answer from below and 
the convergence criterion normally used in the iteration process 
was that the deflections obtained from two successive iterations 
should agree to within a certain limit, in this case five percent. 
When the iteration process is slow, the deflections obtained from two 
successive cycles do not differ very much, but they may still be 
inaccurate compared to the correct final deflection. This is shown in 
Table 6.1 where, for P /P 0.4 and 1/h=20, the error is about 23 per- 
cent. However, this type of inaccuracy will not occur when the process 
converges rapidly. 

In Section 6.4.1 it was suggested that the rate of convergence 
could be used as a check on whether a stability failure is imminent. 
To enable this criterion to be used with Equation (6.31) it is nece- 
ssary to express the convergence check in a different form. If the 
iteration process is considered to have converged when the deflection 
in the i th cycle is within five percent of the final deflection, 


this condition can be expressed by Equation (6.32), 


a 
PA, 
au (6.32) 
I £1.05 ; 
P P i-l 
4, (1+ Wan (7) Ks 
Equation (6.32) may be rearranged to read 
P i Lien 
(7) k, = 0.05 (6.33) 


Thus, the number of iterations required to achieve five percent 


accuracy is given by 
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we 2 0.2e0 (6.34) 


(41) 

log 

ou 

For the particular case where i is not to exceed five, Equation (6.35) 


must be satisfied. 


vu 
> 


ay 
HA 


0755 (6.35) 


a 


Equation (6.35) shows that the deflections obtained in the first- 
order analysis may be used to check on stability failures, and if 
Equation (6.35) is not satisfied, the member stiffnesses may be 
modified until the first-order analysis gives acceptable results 
before proceeding with the calculation of the final second-order 
deflections and stress resultants. 

Alternately, the check on the iteration process may also 
be based on the convergence of the moments rather than the deflections, 


and the convergence criterion may be expressed as 


Q,! + PAs q 

0,1 + PA. 1.05 (6.36) 
1 j 

By substituting for hy and A. from Equations (6.31) and (6.29) in 

Equation (6.36) it reduces to 

i+] | 

PA 2 

= 0.05 (6.37) 

When i = 0, Equation (6.37) is identical to Equation (6.13), and 

thus defines the values of PA,/Q,1 for which a second-order ana- 


lysis may be ommitted. With i = 1 it becomes 
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“9 
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(6.38) 


_ 
A 


0.22 


- 


In this case sufficiently accurate moments will be obtained by 
adding PAs to the first-order moments and the second-order deflect- 
ion bo need not be computed. For values of PA,/Qy1 greater than 
0.22 the second-order deflections must be computed and are easily 


obtained from Equation (6.31). 


6.4.3 The Effect of Sidesway Instability on the Accuracy of The 
PA Analysis 


Methods currently used to predict stability failures 
only consider the bifurcation load which is independent of the 
lateral load. However, when stability failure may occur as a 
result of a reduction in member stiffness, the magnitude of 
the lateral load becomes an important parameter. 

A typical load-deflection diagram for a slender column 
whose load carrying capacity has been sharply reduced due to a red- 
uction in member stiffness is shown in Fig. 6.13(a). This column is 
stable if the lateral load is less than Q1° If the applied lateral 
load is between QV and Q» the design may be unsafe. This will always 
be the case in a PA analysis since the second-order moment will be 
less than My but greater than the actual failure moment. For the 
typical case of a stability failure shown in Fig. 6.5 this would mean 
that the second-order moment lies between the failure curve and the 
short column interaction diagram, and since it is less than My the 
column would normally be considered adequate. 

Fig. 6.13(b) shows a typical load-deflection diagram for 


a material failure. In this case the lateral load carrying capacity 
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is governed by the moment capacity of the cross-section and cannot 
be overestimated unless the lateral stiffness is assumed greater 
than the slope of the line 0-1. The same conclusion can be drawn 
from Fig. 6.3 where the failure curve and the short column inter- 
action diagram coincide. 

An unsafe design does not necessarily arise for stability 
failures in the moment magnifier method, even if the lateral load is 
between Qy and Q, in Fig. 6.13(a), since the moment magnifier is a 
function of Peon: Lf Pes is large enough the magnified moment 
will be greater than My and a column with a larger moment capacity 
is required. This may help in two ways to produce a safe design; 
firstly, the new column may change the failure mode from a stability 
failure to a material failure; secondly, the new lateral load carry- 
ing capacity or may be greater than the applied load. However, there 
is at present no simple method to investigate this type of stability 
failure and the correct failure mode is normally not known to the 
designer. 

If the lateral load carrying capacity could be predicted 
by the use of subassemblage charts, for example, the mode of failure 
would also be known. This will normally involve rather extensive 
computations which makes this method unsuitable for design. This is 
doubly true because the lateral load carrying capacity depends on 
many variables which makes it difficult to establish simple empirical 
expressions or general charts. 

6.4.4 Deflection Limits to Prevent Sway Stability Failures 
Rather than attempting to predict accurately the lateral 


load at which instability occurs, it is simpler to define certain 
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limits which will indicate when a sidesway instability failure is 
likely to occur. Such criteria may be established as a relationship 
between slenderness ratio, axial load and deflection index. 

Fig. 6.14 is a plot of slenderness ratio against the 
deflection index for stability failure for all the columns studied 
in this thesis. It includes data for 217 separate cases involving 
three column cross sections, five slenderness ratios and five differ- 
ent end restraints. The A/L values plotted correspond to the value 
A, in Fig. 6.13(a). This was divided by the real length L of the 
column rather than the effective length kL since sway indices for 
buildings are generally expressed as A/L. Data for columns dev- 
eloping material failure, as shown in Fig. 6.13(b), have not been 
included in this figure. It is seen that the following relationships 
may be used to define values of slenderness ratio and deflection 
index for which stability failures did not occur: 


kL 


A 
([) = Teo0or 


1A 


hs 
fori ee 


90 (6.39) 


HA 


and (A) § 0.005 for kL. 99 (6.40) 
u r 


Equations (6.39) and (6.40) have been plotted in Fig. 6.14. The 
scatter diagram also indicates that for a given slenderness ratio 


there is an upper limit of A/L above which no stability failures 


occurred. However, in these cases the deflections would by far exceed 


acceptable limits and they have therefore not been included. 


Equations (6.39) and (6.40) are based on the effective 


stiffness values computed from the subassemblage analysis. In practice 


the correct effective stiffness is seldom known and thus the effective 
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length factor k is subject to errors. Also, difficulties may arise 
in defining k if it is not obvious whether the frame is free to 
Sway or is braced. 

An alternate approach is to express slenderness by the 
ratio 1/h rather than kl/r. Values of (a/L), for which stability 
failures were not observed may in this case be established from 
Fig. 6.15, which is a plot of ultimate deflections against 1/h for 
217 stability failures and is expressed analytically by Equation 
(6.41). 


Peay 2 Eee for ¢ 2 22 (6.41) 


For 1/h < 22 Equation (6.40) still applies. 

Similar criteria may also be based on service load 
deflections. Figs. 6.16 and 6.17 show plots of minimum and maximum 
service load deflections, respectively, plotted against 1/h for 
the 217 cases of stability failures. These plots show that no stabi- 
lity failures were observed with deflection indices below 1/300 
for minimum service loads and 1/250 for maximum service loads. 

To consider the effect of the axial load, PEE was 
plotted against the smallest value of slenderness ratio and 
deflection index at which a stability failure occurred. These plots 
are shown in Fig. 6.18 and Fig. 6.19, respectively. 

Based on these plots the following empirical expressions 
were chosen to define regions where no stability failures were ob- 


served: 


p 
kL < “u 
ee OW 00s (6.42) 
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(Ay = 0.00095 9.0035 (6.43) 


oc 


P 

u 
oe 0.035 
oO 


These equations have also been plotted in Figs. 6.18 and 6.19. For the 

columns studied, columns having kL/r or (A/L) falling below these 

curves developed material failure rather than stability failure. 
Equations (6.42) and (6.43) may be combined to give a plot of 

Slenderness ratio against deflection index as shown in Fig. 6.20. The 

condition that will have to be satisfied to ensure a material failure may 

then be expressed as 


Ay ¢ __0.00095 


__ 0.00095 ss. 9935 (6.44) 
a Lirgaee 0.02() 


A column with a given kL/r and aden value will be stable if its 

(a/L) is restricted to less than the value given by Equation (6.44). 
Within the range covered by Eqn. (6.44), Equation (6.44) is less conser- 
vative than Equations (6.39) and (6.40). If Equation (6.44) is not satis- 
fied Equations (6.42) and (6.43) may be checked to see if one of them is 
satisfied. 

The results on which Equation (6.44) is based do not include 
values of ails less than 0.1 or greater than 0.8. It is particularly 
important that the equation is not used for P /P<0-1 due to the asympto- 
tic nature of the curve. 

The practical significance of Equations (6.39) to (6.44) is that 
the P-A approximation to the second-order frame analysis can be used for 
frames having A/L values at ultimate which are less than the values given 
by these equations. Based on Fig. 6.14 and Equations (6.39) and (6.40), 
the P-A analysis can be used safely for frames having sway deflection in- 


dices A/L less than 0.005, or 1/200, at ultimate or sway indices less than 
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FIG. 6.18 LOWER VALUES OF SLENDERNESS RATIO FOR STABILITY FAILURES 
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1/300 at service loads. While this conclusion has only been checked 
for three particular column cross-sections in sway frames, the prop- 
erties of sway subassemblage curves suggest that similar limits could 


be used for the design of all practical column cross-sections. 


6.5 Determining the Maximum Moment in Beam-Columns in 
Sway Frames 


The methods of second-order analysis of sway frames dis- 
cussed in previous sections will give the values of the end moments 
acting on the columns. These moments are not always the maximum 
moments in the columns and it is necessary to investigate when the 
end moment is the maximum moment. To do so the moment produced by 
the axial load as a result of the column not being straight must be 
included. This can be done by considering the standard differential 
equation for a beam-column. 

Referring to Fig. 6.21(a) the equation for the moment in 


the column can be written as 


2 
. -_frpdy 
My M. + Py EI 2 (6.45) 


where M. is the moment at the section ignoring the axial load. 


By putting q? = P/EI the equation can be written 


2 M. 
dvy 2 = = I. 
a9 + Qiy I (6.46) 


Differentiating Equation (6.46) twice and substituting dM /dx® for 


(d4y/dx")EI and M, for (d¢y/dx*)EI it takes the form 
2 2 
lon ee 


5+ q MY <7 (6.47) 








dx dx 
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FIG. 6.21 MOMENTS IN A BEAM-COLUMN 
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and the final solution for the moment at any cross-section is 


s ingx +(M,/Mo) sin (ql - a 


M. = My | efral (6.48) 


M) and Mo are the numerically smaller and larger end moments, resp- 
tively. The ratio M1 /M, is negative for double curvature and posi- 
tive for single curvature. 

Fig. 6.21(b) through 6.21(e) shows various shapes of 
the moment diagram that may occur in a beam-column. From these 
diagrams it is clear that if the maximum moment occurs at the end 
of the member the slope of the moment diagram must be positive at 
the end. Thus 


dM. 2 
(--—-) = 0 6.49 
dx iv ts 
x=0 
Differentiating Equation (6.48) and putting x=1 gives 


M 


cos ql 2 a (6.50) 
Z 
for x=0 Equation (6.49) becomes 
< Mp 
cos ql = q— (6.51) 
1 


Equation (6.51) is clearly satisfied for positive Mo/M)- In the case 
of double curvature the equation may be satisfied for Mo/M4 = - 1 but 
cannot be satisfied for Mo/My <-l1. This is true mathematically, 
because the cosine is always greater than minus one and physically 
because the axial load would exceed the critical buckling load. 

Thus, if Equation (6.50) is satisfied the maximum moment will occur 
at one end of the member. 


If the maximum moment occurs between the ends of the beam- 
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column its location is obtained by setting Equation (6.49) equal to 


zero, which yields M 


ae - cos ql 
iy ch 
tan qx = “sin ql (6.52) 


The solution of Equation (6.52) should always be an angle between 
zero and 7, and may be substituted into Equation (6.48) to obtain 
the maximum moment. 

It has been suggested (35) that when the column end moments 
include the PA moments the moment magnifier concept based on a sway 
prevented condition may be used to determine the maximum moment, 

If the moment magnifier turns out to be less than one the maximum 
moment is assumed to occur at the end. 

According to ACI 318-71 the maximum moment is given by 


Equation (6.53), omitting the capacity reduction factor: 


M = M, Ci (6.53) 
P 
u 
L- pt 
cr 
see (6.54) 
where Cn 25096 + 2.4 —= 0.4 ‘ 
2 


In a sway frame, with similar columns which are equally loaded, 

the value of the applied axial load must be less than the critical 
load based on a sway permitted condition and if Ong is evaluated on 
the basis of a sway prevented situation the maximum value of the 


magnifier (for a given Cc.) occurs when 


me) 
sy 


g (6.55 


r 


P 
Pc 


Se) 
TO 


Cc 


where ae and oe are critical loads computed on the basis of sway 
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permitted and sway prevented conditions, respectively. For a given 
column with given end restraints the right hand side of Equation 
(6.55) depends only on the effective length factors, and by using the 
nomographs to determine the effective length factors for various 
values of end restraint it was found that the maximum value 


p S yp D 


cy | Pie can have is about 0.35. If this value is substituted into 


Equation (6.53) it is found that the magnifier is less than one for all 
Mi /M5 less than 0.125, thus indicating that the maximum moment occurs 
at the end for smaller values of Mi/Mo. If the applied axial load is 
only one half of the critical load Peet no magnification would be 
required for M,/M, less than about 0.55. 

A comparison between Equations (6.48) and (6.53) may be 
obtained by assuming values of end restraint and the ratio of applied 
axial load to the actual critical load. Then the parameter ql] may 


be computed as follows: 





qo = (6.56) 
2 
S nm EI 
p?-= (6.57) 
cr (k1)? 
ql=t = (6.58) 
Por 


The range of M,/M, values for which the maximum moment occurs at the 
end is determined from Equation (6.50), and the magnification factor 
is calculated from Equations (6.52) and (6.48) for cases where the 
maximum moment occurs between the ends of the column. 

Figs. 6.22 and 6.23 show how the moment magnifier varies 


with the moment ratio M1 /Mo for various values of end restraint and 
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FIG, 6.22 MOMENT MAGNIFICATION IN BEAM-COLUMNS 
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the ratio oe hd In all cases the restraint has been assumed to be 
the same at both ends of the column. The ratio Pres in Equation 


(6.53) was computed in two different ways: a) P.. =P b 


er cr and b) 


Pen = Pes where Pe is the Euler load. The plots show clearly that 


the use of P_. = P D in Equation (6.52) underestimates the moment 


cr 
magnifier considerably compared to the more accurate solution ob- 
tained from Equation (6.48). This is due to the fact pointed out 
before, that the ratio p /Po will always be low in sway frames. 
A much better approximation is obtained when we = Pes except when 
the column to beam stiffness ratio is very low and p/P a is very high, 
in which case it tends to give results on the high side. A typical 
example of this is shown for peas = 0.9 in Fig. 6.22(a) where 
Equation (6.53) is seen to overestimate the moment magnification by 
a large amount. 

An alternate approximate method may be used to give a better 
estimate of the moment magnifier. First, Equation (6.50) is used to 
determine whether the maximum moment occurs at the end of the column. 
If a magnification is required. the straight line approximation shown 
in Fig. 6.24 may be used. This approach will always give reasonably 
conservative values of the magnification factor. 

The coordinates of points A and B are easily computed for a 
given value of ql. The abscissa of point A is simply, from Equation 
(6.50), equal to cos ql, and the ordinate of point B can be shown 
from Equations (6.52) and (6.48) to be sec(%5). Thus the maximum 


moment may be expressed by Equation (6.59). 
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: 1 
cost 1 
Maye te FRicostaL i cos ql) + 1 (6.59) 


In practice the moment magnifier equation with Sy = Pe 
will yield good results in most cases. Only when the effective 
length factor in the sway permitted condition is close to one, say 
less than about 1.15, will it be preferable to use Equation (6.59). 

The discussion of this section has been based on the assump- 
tion that the frame is free to sway such that the columns bend in 
double curvature. If there is significant lateral restraint provided 
by bracing elements the deflected shape of the columns will approach 
Single curvature. Thus it is appropriate to evaluate ie in Equation 
(6.53) on the basis of the sway prevented nomographs. Since, in 
practice, there is no such thing as a completely braced or unbraced 
frame some judgment must be exercised in each case as to which nomo- 
graph should be used to evaluate the effective length factor. The 
bracing criteria of Sections 6.2 and 6.4.2 may be of some assistance 


in determining the degree of bracing present. 


J 

Tiw © 
4« 
OJ Oh 


sau. at alfds xq 9d at te @L. L duinds nodd 22al 
need 267 not3oee 2 hats to potzzudetb- “ont 
ants jel? dove Yswe oF ort ‘at emett ste “es 


svete! ta sof tine Fe at grorlt +. evuhy 


[oo sdt to sqerde‘b bod af teb ott rman’ 
—_ 


a - ert. 
, ae 
Ths 













f f » . 4 
> ' 7 - >...) = - al 7 a 
Ip 209 =" Tp eo * ae nat” 
: eat } | F, 
Ssups wet tinpsa Sasmom of? sottosng mt - 


i+ notw yfn0 .20269 teom.at eitueay doop bist 


2t notdtbnos bod+tansq (me ont at rotons f 


ot Pe 


_ 


pit 


~ 


fave of stsirqoxage of 3f 2uaT st arn 
i 

sipomon Oss na varg al + its w ef oot 
fae ara 
eieestigits & 26 pnd ‘Aue og bh ova ne a sold 


se60 dowd nf beetowee, ot 2 Swrir: 
| : Tanti th | 
agi vttoatts sad stastane ab tae 2 AGE 


». negara, 


CHAPTER VII 


SUMMARY AND CONCLUSIONS 


7.1 Summary 


In this thesis the second-order analysis of reinforced 
concrete frames has been investigated. 

An analytical model based on the sway subassemblage 
concept was used to study the load deformation response for 690 
combinations of column length, column cross-sections, eight column 
lengths and five beam stiffnesses were included in the investigation. 

The main parameters studied were deflection indices, 
member stiffnesses, stability parameters and second-order moments. 

A computer analysis base on the concepts of column 
deflection curves and sway subassemblages was carried out using 
realistic stress-strain diagrams for steel and concrete and taking 
into account the variation in stiffness along the member caused byy 
cracking. 

The current method of second-order analysis prescribed 
in ACI Code (318-71) and an iterative PA-analysis have been compared 
to the computer analysis to determine their applicability and 
limitations. 

The deflections obtained from the computer analysis were 
used to establish criteria for predicting stability failures. These 
criteria are also the practical limits for the use of the simplified 


design procedures currently in use. 
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7.2 Conclusions 

This study has only considered three column cross-sections 
and only one value of concrete and steel strength. However, since a 
fairly wide range of shapes of M-P-» diagrams was used, the results 
should apply to a wide range of column cross-sections and material 
strengths. 

The analysis showed that the column deflection curve and 
Sway subassemblage concepts may be used to obtain the load- 
deformation response of reinforced concrete frames. However, the 
amount of work involved in this kind of anlysis would in most cases 
make it unsuitable from a practical point of view. 

The deflection analysis of Section 4.5 showed that the 
ultimate deflections corresponding to sway failures are relatively 
large and that in all but a very few cases exceeded current code 
deflection limits. Stability failures experienced the largest || 
deflections and none of the 217 cases occurred at a deflection 
index smaller than 1/200. This indicates that if currently accepted 
deflection limits are satisfied the structure will exhibit a 
material failure. Section 6.4.3 presents some approximate methods 
which can be used to predict the failure mode. 

Service load deflections were also found to exceed values 
specified in most building codes. For high load factors (minimum 
service load) service load deflections varied between one third and 
one half of the ultimate deflection and for small load factors 
(maximum service load) they varied between one half and two thirds. 
Equation (4.7) may be used to estimate the ultimate deflection. 


The moment magnifier and PA methods of analysis were both 
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found to be applicable for material failures, but in the case of 
stability failures they may produce unsafe designs if the value of 
EI is greater than that corresponding to line 0-3 in Fig. 6.13(a). 
In Section 6.4.2 some criteria were established for the use of the 
PA-method, and if they are satisfied the structure will exhibit 

a material failure. By expressing the PA iteration procedure in 

the form of a geometric series it was found that if PA,/Q,1 is less 
than 0.05 the effect of the axial load may Ke ignored altogether 
and if PA)/Q,1 is between 0.05 and 0.22 only the effect of the first 
order deflections need be included. For values greater than 0.22 a 
second-order analysis is required. 

The stiffness analysis of reinforced concrete beams showed 
that the effective stiffness may be assumed to be between 0.5 and 
0.6 times Et: When, in the case of T-sections, the stiffness was 
expressed in terms of the moment of inertia of the web, values of the 


effective stiffness were found to be between 1.0 and 1.2 times El: 
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APPENDIX A 


RESTRAINED COLUMN CURVES 


A.1 Introduction 

The restrained column curves presented in this section 
have been prepared for five different values of the beam stiffness 
parameter K. From left to right the curves are for K = ~, K = 600, 
K = 400, K = 200 and K = 100, respectively. Where there are fewer 
than five curves on the chart, those corresponding to low values 
of beam stiffness are missing, indicating an unstable structure 
for those cases. Thus, if a chart contains three curves they 
correspond to values of K of ~, 600 and 400. For K equal to 200 
and 100 the structure bacame unstable at a very small value of 


lateral load. 
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APPENDIX B 
COMPUTER PROGRAMS 


ed Beam Stiffness Program 


ie a General Description 
This program is written in FORTRAN and 


computes the effective elastic stiffness of reinforced 
concrete beams under gravity load, Lateral load and a 
combination of the two. The variation in stiffness 

along the member is taken into account by dividing the 
beam into a number of short segments which are all of the 
same length. The moment in each segment is taken as the 
mean of the moments at the ends of the segment. The EI 
value of the segment is then determined from the moment- 
curvature diagram of the cross-section. 

When this process has been carried out 
for all segments the beam is loaded by the M/EI diagram 
and the conjugate beam theorem is used to compute the 
flexibility coefficients. The rotational stiffness 
coefficients are then determined by inverting the 
flexibility matrix. 

The beam is assumed to be divided into 
three sections, two end sections and a middle section, 


and the cross-sectional properties are assumed to be the 
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same within any one section. 
When the beam carries gravity load the 
moments produced by these loads are computed assuming 


fixed ends. 


B.1.2 Input Data 
The program can handle any type of 


moment-curvature diagram with a positive slope. However, 
non-linear curves must be approximated by a series of 
straight lines and the coordinates at all points where a 
Change in slope occurs must be read in. The number of 
points must not exceed 99. The coordinates of the origin 


are not read in. 


The load intensity is specified as a 
fraction of the moment capacity of the end segment. For 
gravity loads this fraction is given by the ratio wie/12M, 
where w is the uniformly distributed load, T is the 
length of the beam and Ma is the moment capacity of the 
end segment. For lateral loads it is given by the ratio 
of the sway moment to the ultimate moment of the end 
segment. When the beam is carrying both gravity load 
and lateral load the sum of the two ratios must not exceed 
Des 


All units are in pounds and inches. 
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The following input parameters are required: 


i Type of loading. 
Number of load cases 
Gravity load. 


Lateral load 


ao FSF WW PP 


Number. of segments. All segments 
are of the same length. 
o. Number of segments in left hand 
end section. 
th Number of segments in middle section. 
a. Number of segments in right hand 
end section. 
- Number of points on M-¢ curves. 
10. EI for gross concrete section. 
iia Moment and curvature coordinates 
of M-% curves. 
A detailed description of the input data required 


starts on the next page. 
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Card 
Column 


2-80 


1-2 
3-80 


1.80 


1-5 
6-10 
11-15 
16-20 


21-80 


Name used 
in Program 


CASE 


ALPHA(1I) 


BETA(1) 


NSEG 


NSEG1 


NSEG2 


NSEG3 


DATA DESCRIPTION FORTRAN 
Format 


Type of loading. for gravity I 
load only, 2 for lateral 

load only and 3 for combined 
loading. 


Blank 
Number of load cases. NBETA 12 


Tue 
Blank 


Ratio of fixed end moment SF 10. 


capacity of end segment. 

Omit if CASE=2. Number of 
entries equals NBETA. Combine 
on new card if NBETA >8. 


Ratio of sway moment to Ste 10: 


moment capacity of end 
segment. Number of entries 
equals NBETA. Omit if 
LCASE=1 . Continue on new 
card if NBETA>8. 


Number of segments which the 15 
beam has been divided into. 

Number of segments in left es 
hand end section. 

Number of segments in middle oo 
section. 

Number of segments in right 15 


hand end section. 


Blank 


0 
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Card 
Column 


16-80 


1-10 


11-20 


21-80 


Name used 
in Program 


NP( IJ) 


NP(IJ) 


NP (IJ) 


NP( IJ) 


NP(I,d) 


NP(I,d) 


EM 


PHI 


DATA DESCRIPTION FORTRAN 
Format 
Number of points on M-$ 15 


curve for positive bending 
in left hand end section. 


Same as above; but for I5 
negative bending. 


Blank 


Number of points on M-¢ 15 
curve for positive bending 
in middle section. 


Same as above, but for 15 
negative bending. 


Blank 


Number of points on M-¢ jhe 
curve for positive bending 
in right hand end section. 


Same as above, but for t5 
negative bending. 


Blank 

EI for gross concrete rises 
section. 

Blank 

Moment coordinate of point F10.0 


on M-¢ curve for positive 
bending. Left hand end 
section. 


Curvature coordinate for FLOs0 
same Number of cards= 
number of points on curve. 


Blank 
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B.1.3 Qutput 

The program prints out the stiffness matrix 
and the stiffness coefficients. .It also prints out the eff- 
ective elastic EI value for each end of the beam. The units 


are pounds and inches. 


B.1.4 Listing of Program 
A listing of the computer program starts on the 


next page. 
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DIMENSTON NP(99420)eKETA( 10) ¢EM(9992415) PHI (9962 915)5 

ISMO99 610) ee THETA( 292) 9TM( 99420)4E1 (99920) oS( 292) 5511202) 

GIMENSTON GM(99419)+ALPHAC1O) 

Ce GE nao WEN MES enc ar Ce te CAN LO Nat eh ee ee a a ee 


NOW 





9 C Fee 
190 C €¢*® LCASE=13 GRAVITY LOAO UNLY 
ek: C *#© LCASE=2: LATERAL LOAD ONLY 
12 C #** LCASc=3:2 GRAVITY LOAD AND LATERAL LUAD CUMBINED 
13 C tee 
i — (cs hE he ChS ets 
1s READ (52903) NGETA 
16 IF(LCASE«EO.2) GO TO 1 
ey i 2 C ste 
13 C #** FEAD PATIO CF FeEeMe TO ULTIMATE MOMENT OF END SEGMENT 
19 Cc t* 
S20 RE ATI Se FOS ICRC PRATT oI HT eNGETA) a ee, ae Te 
21 IF(LCAS&eEQel) GO TH 2 
22 C +t 
23 C #8¥ READ PATIO OF SWAY MOMENT TO ULTIMATE MOMENT OF CND SEGMENT az 
246 C tee 
25 C #$# NOTE THAT ALPHA(L)+BETAC(L) MUST BE LESS OR EQUAL TO 1.20 
Ss eee 
27 1 READ( 56904) CHETAC LE) o L=) eNLTA)D 
28 2 READ(5+900) NSCGeNSEGI »sNSEG2,NSEG3 
a. -’ C ### 
34 C #¢*% FREAD NOs GF POINTS GN M=PHI CURVES. NP(TeJ)de 
31 C ##® [ IS THE SEGMENT NUMBER» - 
~~ sa, -— G& ** Jer RERENS TO THE M=PHT CUPVE FUR POSITIVE SCENOTNG,  — > 
33 C *## J=2 REFERS TO THE M-PHI CURVE FOR NEGATIVE SENDING. 
34 C e*e 
is = eGONTAO Sedietedaai = 
36 READ(5e90C) (NP(155) eJ=l 42) 
37 RLAD656499) CNPENSEGI 41 0d) ey=le2) 
Si 2 mum am 2 aac a aaa 
39 READ(S5+906) (NP(NsJ) eJ=1 02) 
43 059 39 T=2.NSEGL 
—. ee CFs PANE) eh i ° P : or aad 
42 39 MPC 167 d=NP(122) 
43 N=NSEGI+2 
SS NEN SEC NSE eee 
495 00 40 [=N.NN 
46 NP CT +1) =NPCN—151) 
eT —~ 40 NP(1Ts2)=NP(N=-1 92) : sea 2 
43 NNN=NN4#2 
49 DO 41 L=NNNeNSEG 
——— eS NCL OE SSE NTA pr) 
si 41 NPC1s2)=NPCNN41 52) 
s2 READ(S5+913) ECIG 
=1S53°-- ¢ ##% na 
s¢ C @#¢# PEAD M=PHI COORDINATES (MOMENTS IN LG—IN) 
$s Cc ¥st 
Se CaS ta ee 
s7 NPL=NP(1 04) 
Ss DO 3 K=1.NPL 
ss ‘ 3B READ(5.904) EM( LedeK)ePHIC1e Sek) 
60 00 42 I[=2-NSEGI 


61 DO 42 J=1s2 
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62 NPL=NP(I.J) 
o3 DO 42 K=1eNP1 
64 EMC Le de K VEEM( Lede K) 
65 42 PHI (Led eK)=PHI (Le deK) 
66 DO 43 J=1e2 
o7 NPI=NP (NSEGI41.4N) 
ee _c ihel. 
69 43 READ(5.904) EM(NSEG1L 41eJeK) ePHI(NSEGI 41 eJeK) f 
7e O00 @& [=N.«NN 
ee) gee DO 44 J=1.2 
72 NPL=NP (IJ) 
?3 OC 44 K=1,.NP1 
eT TICE EM 5 ACLS ee ae it eS, a ae DE Lie a es Fes 
75 44 PHI (1.5 eK)=PHICNSEGL F1e30K) 
76 00 45 J=1.2 
ys aia NPL=NP (NSEG-NSEG34164) > a 
73 DO 45 K=1,NP1 
79 @5 PEAD(5+998) EMI NSEG-NSEG341 040K) o PHI (NSEG-NSEG 341 69K) 
ee NEN SLE LENSe Gets. 
gL DO 46 IL=N2.NSEG 
382 20 40 J=142 
+35 "= ~~ NPLENPC(T su) Se te <a ; : es Soo ee eee 
cry DO 46 K=1.sNP1 
85 EM( 1. JsK)=EM(NSEGIFNSEG2Z F165 0K) 
—S6 ES PRIUS eX TEPATUNSEGT INSEGZ FIs eK) 
RT Cc eee 
as C @*® READ LENGTH OF BEAM (8L) IN INCHES 
3g ——~"-C ee z sid 
29 READ (5.902) AL 
91 WRITE(0+921) 
Sa ee eR TEGO .oeek en ere oe 
93 WRITE (64905) LCASE 
94 IF(LCASE*EQ-2) GO TU 31 : 
3s WRITI(6.927) : 7 > ae ee oe Me eg Gaadae ee Ee eet ate 
96 WRITE(6.907) (ALPHA( I)» I=1eNJETA) 
9? IF(LCASE.sEQe1) GO TO 33 


9? WRETE(6,9C07) (SETACI) sl=1 eNBETA) 


190 38 wRITT(6.923)9L 
101 WPlLTE(6e903) NSEG a = = = = ae 
1c2 WRITE(G+.914)ECIG 
103 WRITE(GO+912) 
‘ —_ - -* 
10s WRITE(6 0925) 
106 WRITE(6-911) 
107 3 EGL=BL/NSEG 
1o3a X=06 
109 IF(LCASE*«EQe2) GO TO 4 
ee ee Or ee Se OE AS ee er eae fey ee ee 
ili X=0.0 . 
112 NPIL=NP(1.2) 
113 GL=12.C FALPHA(K) FEM(1.62.NP1L)/( SL **2) 
114 GML=-GL#®OL*#2/12.6 
115 X=X*#SEGL 
Sy ene seeped SR USS i, Lok ee 7 ae aaa a a a a a a a a a a a a ee Tae Tea 
117 GM2=GL*(6 00 BL EX—CL ##2— 60 0 EXH#2)/1260 
113 GM(I4K)=(GM1+#GM2)/220 
m9 X=X+SEGUL 
120 S GM1=GM2 


i2i OO 81 L=leNSETA 


tr 
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122 Bl WRITE(Gs931) GAC 1 +L) + GM(NSEGoL) 
123 4 IF(LCASE.EQe1) GO TO 6 
124 K=NP (12) 
125 OO 7 L=1+eNBETA 
126 X=SEGL 
127 SMI=CBETACL)*EM(1+25K) 
ito os Cae ee hh a me Sa... =. +  — = ee 
129 SM2=HETA(L) #OM( 1 620K) #4 (260*(BL—-X )/BL)—-160) 
130 SMC IT eLI=(SM1L4+SM2)72.20 
~ ASI _— SM1=SM2 
132 7 X=X+SEGL 
133 IF(LCASE+GTel, GO TO 9 
i SR OR CDN Dc ng ee ee 
13s 00 19 L=1lseNSEG 
136 10 TM(IsL)J=GM(IeL) 
=137F Tria GO TO 14 = 
138 9 IF (LCASE-EQ-3) GO TO 11 
139 DO 12 I=1+NSEG 
Le STR g pS age 21 eae oe ay at FT ig a agg arms a eg ei ce a mA a 
141 12 TM(IsJI=SM( Ie) 
142 GO TO 14 
~18a3- — It DO 13 I=1eNSEG 
144 OO 13 J=1,NNETA 
145 13 TMC EL 2 JV =GM( I+ J) *+SM(1 9d) 
—1456-—-"—cC” ¥e8 scale ee eit —s SS Ser amnion £5) 
147 C @** COMPUTE EI FOR EACH ELEMENT 
143 C ste 
~149 ~~ 61 M=0 
iso 14 IF (LCASE.GTe1) GO TU 15 
151 DO 16 L=leNJETA 
LL °° °«&° « °« °«>I°xD LOD LS NDEGo Fi FO GE TREE le RR a a aA eae ee Tee ec ee See 
193 J=2 
154 TFC TMC 1 eh) eGT 2005-0) J=1 
135 IF CABSC TMC Loh D2 eGT EMCEE e010) GO TO 17 . ; : ee =. oe, 
156 NL=NP( 14 J) 
197 IF CC TMC LT eL IP -EM(C Le SeNL D2 eGT*1020) GO TO 200 
ee A ee 
1s9 ECC L sL=EMC LT eJe lL DSPHI CE sJol) 
160 GO TO 16 . 
161 — 200 EICIsLI=90.-N09NGNGOOOCONOL > 7 a a are SA ee Oia oR C 
162 WRITE (62933) 1,B8ETA(L) eALPHA(L) 
163 GC TO 16 
ee ee en eee eee ee ee ek eo ee oe ee 
165 19 TFCABS(TM( Isl )) eGToEM(ITe4eN))? GO TO 18 
166 ELC sLV=C EMC Lede NI-EM( TL ep deNm1L)IS(CPHIC Lede NI-PHI( Lede N-1)) 
167 GE Ee 
108 18 N=Nel1 
169 IF(NeGTeN1) GO TO 203 
ee ——————————————eeeeeeeeeeeee 
7. 203 EIC(IsL)=0-000000000001 
172 16 CONTINUE 
173 == = GO TO 34 - : > . = ee 
174 15 DO 21 L=1eNGETA 
“7S OO 21 I=1+NSEG 
=: aT - 
177 TFC TMC 1 oh 2 eGTeO0e) J=1 
178 NI=NP(I.J) 
179 de TF CO TMC IL se LP-EMOC1L ep deN1 D0) eGT 210-60) GO TO 201 ES ene oe 
180 M=L° 


161 TFC ABSCTM( IT oh )) eGTeEM(I sJe1)) GO TO 22 
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182 ELCIieLY=EMC Te Joi D/PHI(IT eset?) 
183 GO TO 21 
184 201 ELC TeLI=0.00000000000001 
185 WRITE (66933) 1 sSETACL F,ALPHA(L) 
180 Go TO 21 
187 22 N=2 
—_ aa. — 5» a2, Ot ROSti Matshita cleat lel iti) 60 10 2s°”6~CCSC™C~*C:C 
189 ELC eLY=CEMC Lede NI—-EM(C IE pS eN-LIISCPHI (Led eNI-PHIC Le JeN-1)) x 
190 GO TO 21 
“191 ——""" "24 N=NF1 7 ital % ; : a ee — sain 
i192 LF (NeSGTeN1) GO TO 202 
193 GO TO 23 
? =0e0 ayes Gee Ge a ee ee ee ae ee 
195 21 CONTINUE 
196 34 CONTINUE 
197 TF (KOUNT.2EO.0) GO TO 72 ‘ s ai ee aoe 
198 CO 71 L=1isNSETA 
199 DO 71i 1=1.NSEG 
ei. Sein tw Ce CREM CE SET) et tte ei Cl pik 
201 71 CONTINUES 
202 72 CONTINUE i 
203 3 KOUNT=KOUNT#1 : = a a ee ee Sa a a a ns 
204 CALL ANCOL(SEGL eE Le NBETAsNSEG»s TMeo lL CASE o Gl oe SMe GMeBL sALPHAs NP 
26s 1 EM) : 
a Re Nanna TET ED aac ROR Od EE ee ee ee Se Se 
207 DO 70 L=1.NSCG 
203 TO SPC Qe deb Ces} 
Soe IF(KOUNT.LT.5} GO TO ot i : ; ‘ y i. ae 
Z19 Cc eee 
2:1 C =** COMPUTE OVERALL BEAM STIFFNESS USING CONJUGATE HEAM METHOD 
a ee eee eS ee ea core ee nie 
213 Cc #8 
21la C ¢#* UNIT MOMENT AT LEFT END (1) . 
i150 CCCU SRE ia ; F a rie ies eR Se ee ee a Pe oe ee ee 
2:6 OC 29 L=1eM 
217 THE TA(121)=920 
hee Oe ete Tk eee tee eee 
219 THE TA(2.2)=020 
220 THE TAC 122)=0-0 
221 X=SCGL/2.0 
222 OO 26 I=1sNSFEG 
223 THETACL Ss LI=THETAC 1 ol 4( ((BL-X)/JBL) MH#2)e(SEGLZEL( 1 +L) ) 
ee EMC PRL SU ED STL At eek Po 88 COL RIF OLS PL AZOLE SC OLoCLsEcLtCL.sclF.6C™C™C™C™C™C™CCO,YSCCNC#*“‘##§$NSNNSNNNNNN"-S 
225s 26 X=X+*+SEGL 
226 Cc #2 
227 C #8 UNIT MOMENT AT RIGHT END (2) 
223 C ee 
229 K=SEGL/2-0 
eee eee EIEe Hee eee Se eas ee a ee ee Ne ee ee ee 
231 THETA( 202) =THETAL2 62) 4( (X/3L) #2) ®(SEGLZEI(1.L)) 
AEF 4 THE TAC 1 s2) =THETACL 2 )—€ (X4BL ) #( BL =X) /BLIS(SEGL/EI(I4L)) 
233 27 X=X+SEGL 
234 DETHATHETACL se LI *®THETAC 202 )—-THETA( S61) ®THETA( 102) 
235 SC(telLY=HTHETAC2].2)/0ET 
ere EE eae ee A pI ee SS 
237 S(2.1L)J=-THETAC2<1)/0ET 
238 $(2e2)=THETA( 1.2 )JOET 
239 OO 28 L=te2 
240 OO 28 J=1-2 


261 23 SICL6JI=ASTI.JIOGL/SECIG 
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242 WPITE (6-925) 
243 IF(LCASE«EQe1) GU TO 37 
244 TF(LCASE*EO63) GG TO 36 
245 WRITE(6+¢918)4ETA(L) 
246 GO TO 30 
247 36 RA =ALPHA(L)/HFTA(L) 
_ a. . _ anIR Ete dl SORA Le Eee PO ee 
249 GO TO 30 
250 37 WRITE (6.920) ALIVHA(L) 
~ 2am "30 WRITE(6s9039) (CSCI ed) eJ=ls2)eT=162) 
252 WRITE( 65915) 
253 


254 @RITECG.I2ZT 


WRITE(66916) ‘(SIC Ted) oe J=1ee) ol=1e2) 








255 SLOV4=S(1.1)¥*3L/74.0 
256 WRITE (Ge923)SLOV4 
257 + ay SLOV4=S(242)*BL/420 
256 29 wRITE(6+939) SLOVS 

259 DO SO L=1sN3ETA 


zeg WRITETOSII2ZD 


261 @RT TO (62907) GETA(L) 

262 WRITE (CE +931) (GM(IsL) sT=1sNSEG) 

263 WRITE(4£5932) y oo yeee a a i 
264 ALTE (64¢931) CELC Teh) sL=1leNSEG) 

265 50 CONTINUE 


rn ie eek ee ee eae ee os ee 


267 20 WRITECG¢917)I+BETA(L oJ 
268 IF(MeGTel) GO TO 34 
~ 269 ~~~ 900 FORMAT(I6I15) : gc RS ea 
270 901 FGRMAT(I1) 
27) 902 FGORMAT(FI0O.0) s 
ase = SOs EORMAT UL ET re ae eee 
273 904 FORMAT(8F10.0) 
274 905 FORMAT(IHO.SXe»* TYPE JF LOADING =*,12) 
"27S 5 GOT FORMAT (1HO65X 1 0FS3e2) . Nes Tae 
276 908 FORMAT(1LHOsS5Xs*HKEAM I[S DIVIDED INTO *4I136* SEGMENTS!*) 
277 909 FORMAT(1HO 5X pl PZ2EI1S 24) 


wee SP near Oldie ae Loe MOM nhs. SSSSSC*“#EN-—— Se 
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91l FCRMATCIHOsSXs*ROTATIGNAL STIFFNESS COcFFICIENTS: *) 
260 912 FOPMAT(1IH1) 
“281 ~~ 913 FORMAT(F1IS5.-5) <7 ie © eee 
262 F914 FOPMATCIHOe DXe*ECIG=" oF lSelse* LU-INe#2*) 
2383 91S FOSMATC1IHO+SXe* STIFFNESS COEFFICIENTS EXPRESSED AS K=SL/EI1:*) 
—2355 —_ST6_ FGRMATCUTIAO SX eo IP2CISe4) : Oi aaa a 
285 917 FORMAT(IH1s5Xe*THE MOMENT IN SEGMENT NOe* eI 30% 1S GREATER THEN THE 
286 AULTIMATE MOMENT OF THE SEGMENTs BETA=*eF4e2s*GBENDING CODE="4I12) 
267 ——~ GIB FORMATCIHOsSXse* BETA =",FSe2) 
238 919 FORMATCLHO+SXo*F eEeMe/SWAY MOMENT =" oF S02) 
289 920 FOPMAT(IHOs5SXs*ALPHA =" +FSe2) 
—230 2 t-FTOrPtatTtinorsxe* tor Poe 
291 922 FOPMAT (6X%0° -------—-—* ) 
292 923 FORMAT(LHO+SXe "THE LENGTH OF THIS BEAM IS ".FSele* IN*) 
293 . 924 FOPMAT(SXs*O UT PU T*) 
294 925 FORMAT (SX +s %=——<-—--= ee SD 
295 926 FORPMAT(1HO) 
i a Fi aad Ee phy Ee Mh er AE mh Ug DAG A 1 ai aa ac a 
297 928 FORMATCIHO+1OXs+*LEFT END' el PELS.S) 
298 929 FORMAT(1HO»sSXe*EFFECTIVE ELASTIC EI VALUES:*) 
299 ; 930 FORMAT(IHOI,LOKs* RIGHT END* s+ IPELD6S) 
300 931 FORMATCLHOs 1PGE10 22) 


301 932 FORMAT(/S///) 
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3R5 tole A aint Pha Te MD aL tl A ELEMENT® -134* IS ZERO FOR GETA =", 
304 stop 
305 ENO 
306 SUBROUTINE ANCOL (CSEGL sE feNGETA+NSEGs TMeL CASE eGlhe SMe Ge BL sALPHAsNP 
307 1e&M) 
STSCI ENS TON ACPAATIUTINPUIS SOV eEM( 9S Sold) 
309 DIMENSION SM(99 410) + TM( 99920) sE1 (99420) »GM( 99010) 
310 DIMENSION FEM(5450) 
io cy ce DIMENSION AO(50)+XBAR(50)+P(50) eMY(50)+TY(50) 
312 PEAL MYelY 
ais € ®**® COMPUTE CENTROID OF ANALOGOUS COLUMNe AREA ANDO LOAD 
ee Gye. Gets Noe? A tac tl 
315 X=SEGL/240 
316 A=0-0 
317 B=0.90 ee nate i 4 3 
313 C=0-0 
319 NPL=NP (1.2) 
eevee a ot fate AER MCL SeeNrLio koh eta. ek ee) eh ee 
321 OO 1 IT=il+eNSEG 
322 A=A4SEGLZEI( 1 ol) ; 
Fe Moses mete "S234 SEGL/EF CE 0k) ~ 4 ; . ea) pe eee 3 - ora a - 
3e3 C=C+# (GL *X*(BL—-X) *SEGL/S( 200 EI (LoL) )) 
325 1 X=X+4#SIGL 
sats. 7... SDAA 
327 AO(LI=A 
328 P(Ld=C 
329 3 WRITECE6, 102) XSAR(L)+AO(L )-P(L) 
3390 TO] FURMATC1H00S5X41P3F1063) 
331 2 CONTINUE 
eee COMPU Te Tt Ao 
333 DO 3 L=Il-eNNETA 
334 AS=0.0 ; 
as" &8=0.0 i 3 ae ; = eye sg i a haa a oo a eee ey te 
336 X=SEGL 
337 NPL=NP(1-2) 
Ata a i. mie © Cashes SACPNAC EL. SEM Lycs Ne Liv tOoLssel>..c.- °° #-- - 
339 DO 4 I=leNSEG 
340 AA=AAF (GL*X* (BL -K)72 260) ©(CX-XBAR(ILIIZET (I «tL? 
331 BS=BS*+SEGL#*#3.0/( CIC 1 +LI*#12.0) + SEGL*¥(X-XBAR(L) )*#*2/ET (1 2L) 
342 4 X=X+SFGL 
343 MY(L)Y=AA 
a ee er ee NOL, Bakeware ae a Ree ee eee Se arr) ie oie ee SM ee ee 
345 3 CONTINUE 
346 C #@* COMOUTE FIXED END MOMENTS 
347 00 5S L=1.,NHETA 
348 FEM(1LeL)=PCLISAOCLI FMY(L) ®(—-XBAR(LI)DJIY(L) 
349 S FEeM(2eL)=P(LISAO(L) #MY(L) #(GL—XBAR(L))ZIYVIL) 
i ee ee ee eee Sr ek Se ne 6 ae. a ae ee it — 
351 6 WRITE (G*e1 00) FEMC Leh) eFEM(26L) 
352 TF(LCASE-EOe2) GO TO 12 
353 OC 7 L=1eNJETA 
353 GM1=0.6 
355 X=SFGL 
Sap a ah ea Pe PL Se ta IDM Ge SL Sip a et a ian Ama ee Rg aa cama NR a Na a a RR 
357 AGQ=AESCFEM(2sL)) 
353 NPL=NP (1-2) 
359. « GL=12-.C#ALPHA(L)I*EM(1.2eNPL)/ (GOL 8 #2) 
360 DO 3 T=1L«NSEG 


3ol GMA =SLEKECIL =X) 7260 
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302 GM( ToL) =(CM14+GM2 97220-( ASL + (AB2—-ABL ) ®(X-SEGLZ2-0 )ZAL) 
363 X=X+*#SEGL 
364 8 GM1=GM2 
365 7 CONTINUE 
366 TF (LCASE.EQ.23) GO TO 10 
307 BDO 9 L=leNSETA 
= = soar. OO ae bel .NseG 
309 9 IMCI+LI=GM(I6L) 
370 GO TO 12 
371 10 DO Il L=1leNSETA 
372 OO 11 I=lsNseG 
373 Ll TMCT eLI=SMOCLeLIFGM(IeL) 
_a34 Ohba cj LeccumT Inve ; eer 
CW a 101 FORMAT CIHO SSX s1P2E1 063) 
376 RETURN 
377 ENO 
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B.2 Restrained Column Curves and CDC Program 


B.2.1 General Description 
The program is written in the FORTRAN 


programming language. 

The first part of the program generates 
a set of Column Deflection Curves for a given axial load 
and all the geometrical data of the curves are stored. 

In the second part of the program the data 
of the Column Deflection Curves are used to compute restrain- 
ed column curves for five values of beam restraint. 

When the first two parts have been carried 
out for all values of axial load specified, the program 
goes on to compute effective elastic EI values, effective 
length factors and the ratio of the total end moment and the 


PA-moment. 


B.2.2 Input Data 
The program can handle any type of M-P-9 


diagram with a positive slope. However, non-linear curves 
must be approximated by a series of straight lines and the 
coordinates of all points where a change in slope occurs 
must be read in. The last point on M-P-% curves terminating 
with a negative slope is taken as the point of zero slope. 
The number of such points must not exceed 24 per curve. The 
coordinates of the origin are not read in. 

Two service load conditions may be analysed 
with each ultimate load case. The load ratios are read in 


using an implied DO-loop and the sequence must be as follows: 
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ultimate 
load case ] maximum service load 


minimum service load 


ultimate 
load case 2 maximum service load 


minimum service load 


and so on. Thus,fields 1,4,7 etc. are reserved for ult- 

imate loads, 2,5,8 etc. are reserved for maximum service 

load and minimum service load are entered in fields 3,6,9 
etc. The subscripts for ultimate load are then (1+3n), for 
maximum service loads the subscripts are (2+3n) and for min- 
imum service load (3+3n) where n is an integer ranging from 
One to the number of load cases entered. If, for a particular 
load case, values for one or two of the three load ratios are 
not to be entered the appropriate fields should be left blank 
and the input parameters START,END,and STEP should be given 
values such that the program does not use these load ratios 
during execution. START, END and STEP specify the subscripts 
of the load ratios to be used during execution of the prog- 
ram. If eight load cases are to be investigated there can be 
as many as 24 load ratios. To get the results for ultimate 
load for load cases one through four START must be equal to 
one, END must be equal to (1+3x4)-3=10, and STEP equals 
three. The program will then use the load ratios given the 
subscripts 1,4,7 and 10. If only mimimum service loads were 
to be investigated the values would be 3,(3+3x4)-3=12 and 3, 


respectively. 
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If the load factors are different from 1.70 


and 1.275 lines 115 and 113 must be changed accordingly. 


Listed below is- an outline of the input 


data required. 


Li 


o0oclUCU UOOlmUlUlCOUCUCUMNC CDN Cie 


Parameters to specify for which axial loads the 
computations are to be carried out. 

CDC parameters: initial angle, angle increment 
and segment length. 

Column properties: type of column, material prop- 
erties, steel percentage and cover. 

Number of M-P-# curves. 

Number of points on M-P-@ curves. 


Axial load ratios. 


- Moment coordinates of all M-P-¢o curves. 


Curvature coordinates of all M-P-# curves. 
Gross area of concrete cross-section. 


Values of beam restraints. 


A detailed description of the input data 


required begins on the next page. 
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DATA DESCRIPTION 


The number of the load 
ratio at which the comp- 
utations are to start. 
The number of the load 
ratio at which the comp- 
utations are to end. 
Determines which load 
ratios are to be excluded 
from the computations. 
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Initial angle of CDC, radians 
Blank 
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Compressive strength of 
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Yield strength of steel, psi. 
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depth. 
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Column 


1-80 
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1-80 


1-80 


2-9 


10-80 


2-51 


Name used 
in program 


NP] 


PRAT 


EM 


PHI 


AG 


CAY 


Bizwa UNtpUut 


The program prints out the coordinates of the 


DATA DESCRIPTION 


Number of points on 
M-P-¢ curves. Continue 
On next card if N@>16. 


Values of axial load 
ratios. Continue on next 
card if N@>8. 


Moment coordinates of 
M-P-» curves. Continue on 
next card if NP1>8. Repeat 
for all curves. 


Curvature coordinates of 
M-P-o curves. Continue on 
next card if NP1>8. Repeat 
for all curves. 

Blank 


Gross area of concrete 
cross-section. 


Blank 


Blank 


Beam restraints 


FORTRAN 
Format 


1615 


8F10.3 


8F10.3 


Se 0.3 


FE.5 


SLO sO 


restrained column curves for ultimate loads on output unit 


6. Service load parameters are printed on output unit 2. 


Effective EI values, effective length factors, deflection 


index, ultimate lateral load etc. are printed on output 


untt 3. 


The coordinates of the column deflection curves 
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are not printed out. 


B.2.4 Listing of Program 


A listing of the computer program starts on the 


next page. 
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